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DES MOLÉCULES POLYATOMIQUES 
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Chaire de Physique de l'École Supérieure de l'Agronomie, Szczecin 


(Reçu le 14 Juillet 1959) 


On sait [1] que la solution formelle du probléme de l'anharmonicité des molé- 
cules polyatomiques consiste à la présentation du terme d'énergie sous la forme de 
série, où on rencontre les fréquences harmoniques w; et les constantes anharmoniques 
x, (en négligeant les termes suivants). En principe, elles sont diverses pour des 
molécules isotopiques de la méme combinaison. Comme leur désignation exige la 
connaissance de beaucoup de fréquences du spectre de vibration, les fréquences 
harmoniques sont déterminées avec assez d'assurance pour un petit nombre de molé- 
cules triatomiques. 

A cause de la difficulté du calcul des fréquences harmoniques des molécules 
polyatomiques, divers auteurs emploient des méthodes approximatives, qui se rédui 
sent à deux sortes: 

a) introduction de ,la masse spectroscopique^ des atomes et ; 
b) emploi de diverses formules, qu’on peut amener à la formule de la fréquence 
fondamentale de la molécule diatomique 


v = w (14-50) (1) 


Les méthodes de la première sorte sont des méthodes pas trop exactes et peuvent 
servir en premier lieu pour tenir compte de plus grandes anharmonicités liées avec 
les vibrations de valence C—H et C—D (bien que Maiantz [2] a montré que ces 
méthodes usées conformément sont de méme importance que les méthodes de la 
deuxième sorte; mais à cause des calcules trop compliqués elles ne sont pas commo- 
des). Les méthodes du deuxième genre, employées conséquemment pour toutes 
les vibrations exigent l'introduction de trop de constantes d'anharmonicité et peuvent 
étre appliquées seulement pour les molécules de la méme symétrie. 

Maiantz [2] a proposé aussi la formule (1), mais il a fait une certaine essentielle 
supposition supplémentaire. La constante d'anharmonicité n est, chez lui, commune 
pas seulement pour les fréquences des molécules isotopiques de méme symétrie, 
mais pour toutes les fréquences, qui sont attribuées en sens approximatif à des mémes 
coordonnées naturelles sans égard pour la symétrie d'oscillation et pour la symétrie 
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de la molécule isotopique. Le nombre des constantes d'anharmonicité n'est pas alors 
grand et leur addition comme inconnues ne complique pas le calcul des constantes 
de force fait à base des fréquences expérimentales de molécules isotopiques. P. e. 
Maiantz dans son travail [2] désigne 5 constantes de force et 2 constantes d'anhar- 
monicité (7, pour oscillations de valence C—H et C—D et 7, pour oscillations de 
déformation) de la molécule du méthane et de ses dérivés deuterées à base de 26 
fréquences en appliquant la méthode des plus petits carrés. C'est dommage qu'on 
ne connait pas les fréquences harmoniques du méthane et à cause de celà on ne peut 
faire une assez certaine conclusion sur la possibilité de la calculation des fréquences 
harmoniques. 

L'auteur de cette notice croyait qu'il était indispensable de faire un calcul ana- 
logique d'une molécule avec les fréquences harmoniques connues de l'analyse d'un 
vaste spectre d'oscillation. | 

Les résultats du calcul des constantes de force et des constantes d'anharmonicité 
fait pour l'eau et ses dérivés deutérées aprés la supposition de Maiantz par la méthode 
des plus petits carrés à base des fréquences fondamentales présentées par Benedict, 
Gailar et Plyler [4] sont placés dans les tableaux 1 et 2. Dans le tableau 1 nous voyons 
les constantes de force et d'anharmonicité avec des erreurs statistiques et par com- 


Tableau 1 


Constante Calcul actuel Aprés Darling et Dennison [7] 


K, (8,441+ 0,022) - 105 dynes/cm 8,428 - 105 dynes/cm 
K, (0,70280,0056) - 105 erg/rad? 0,7044, - 105 erg/rad? 
h — (0,0981 + 0,0033) - 105 dynes/cm —0,1051 : 105 dynes/cm 
a (0,2422+ 0,0434) - 105 dynes/rad 0,2415, - 105 dynes/rad 
UP — (0,011834- 0,00032) - 107? cm 

Ny —(0,02120+0,00277) - 10-3 cm 


Tableau 2 


Calcul acmel Après Darling Après Benedict, Après 
et Dennison [7] | Gailar et Plyler [4] | Khatchkuruzow [6] 


3829,9+5,1 3825,32 ` 3832,17 3835,37 
1651,4+8,1 1653,91 1648,47 \ 1647,59 
3939,1+5,4 3935,59 3942,53 3938,74 


2761,4+2,6 2758,06 2783,80* 
1208,3-- 4,2 1210,25 1206,39. 
2886,2+ 2,8 2883,79 2888,78 


2821,7+2,6 = 2824,32 
1447,2+6,2 i 1440,21 
3886,8+5,3 3889,84 
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paraison les constantes de force de Darling et Dennison [7] (les symboles aprés Elia- 
chévitch [3]). Dans le tableau.2 on voit les comparaisons des fréquences harmoniques 
calculées (avec des erreurs statistiques) et des fréquences harmoniques obtenues 
par d'autres auteurs à base d'analyse du spectre d'oscillation. Le seul plus grand 
désaccord, dans le lieu marqué d'étoile, est, évidemment, une faute d'imprimerie 
dans le travail [4] (les auteurs écrivent, que leurs fréquences réalisent la régle du 
produit parfaitement, mais la substitution des fréquences montre un désaccord 
de 2%). On a calculé aussi les fréquences harmoniques et prévu les fréquences fon- . 
damentales pour molécules avec tritium en obtenant un résultat conforme aux 
expériences de Staats, Morgan et Goldstein [5] (tableau 3). 


Tableau 3 


7 Fréquences Fréquences 
Fréquences 


à fondamentales fondamentales 
harmoniques 


calculées observées 


2367,1 2300,8-- 1,8 
1372,5 1332,6-- 5,2 
3885,7 3707,14-4,8 


2359,4 2293,5-- 1,8 
1116,8 1090,3+ 3,5 
2832,7 2737,8+2,5 


2299,3 2236;7+ 1,7 
1016,8 994,9+ 2,9 995,9 
2437,1 2366,8+1,9 2370+5 


Il est difficile d'admettre le calcul ci-joint, malgré ses résultats, comme preuve 
de raison de la supposition de Maïantz. C'est seulement la vérification faite sur un 
exemple; son succès fait cette supposition vraisemblable à certain degré. C’est impor- 
tant à cause de l’éventuelle possibilité de désigner plus exactement les constantes 
de force des molécules polyatomiques et d’apprécier leur fréquences harmoniques, 
ce qui d'autre côté peut. faciliter linterprétation des fréquences composées. 

L'auteur désire d'exprimez sa reconaissance à M. Thomas Daca pour des nom- 


breuses et précieuses discussions. 
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An analysis of 53 40 — p + rn decays found in a stack of emulsion exposed to 
a beam of K- mesons gives Q4 = (37.58+0.18) MeV and M4 = (1115.42:£0.19) MeV. 
From an analysis of the results of other laboratories on the standard range of protons 
from the decay of 2344 hyperons at rest and the results obtained in the calibration of the 
emulsion used in the present work one obtains as the best estimates the value R,, = 


= (1678.6+3.2) u, which corresponds to My, = (1189.43:- 0.31) MeV. 


I. Experimental pari and scanning 


` .The energy of the decay 4? — p--z- has been measured in recent years by 
several groups of investigators using chamber and emulsion technique. Results 
obtained by the various groups show discrepancies which in some cases are large 
in comparison with the error quoted. Using nuclear emulsion, one might expect 
to get in this case an accuracy somewhat better than 0.2 MeV. Such an accuracy 
is of interest, e. g. in evaluating the binding energies of the 4° hyperons in hyper- 
fragments, especially in the lightest ones. It was thus felt that further efforts in deter- 
mining the mass of the 4° based on larger statistics would be justified. 

As a source of AP hyperons we chose interactions of K- mesons in nuclear emul- 
sion. We used in this work a stack of 180 plates 10x 10 cm x 6004 of Ilford G5 emul- 
sion exposed to the enriched K- beam (~ 300 MeV/c) from the Berkeley bevatron. 
We expected to find there a rather large number of A? decays with a relatively low 
background of pseudo A9 — p+r events. 

Eighty pellicles from the central part of the stack were available to us for this 
analysis. A strip 2.5 cm wide in which most of the kaons were expected to stop was 
area-scanned for o,- and og- stars. An average of about 62 oy. stars and about 
105 c,. stars per plate were found in the scanned area. An additional scanning for 
purposes of calibration provided us with 22 2 decays at rest, ; 
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In order to find the 4? decays the tracks of light mesons ending in a ø star were 
followed back for a distance up to about 20 mm. Of the 680 tracks traced back in 
this way 54.5 % were abandoned after about 20 mm, 2.7% originated from stars 
with one prong, 1.9 %were due to 7” mesons emitted from oc stars, 11.3 % originated 
from stars with three or more prongs, 10.4 96 traversed the emulsion prior to or 
after exposure of the stack (no prolongation in the neighbouring plate), 10.5 % were 
abandoned before the range of 20 mm was reached, as they originated outside our 
part of the stack, came from a region of high distortion, or were lost during the pro- 
longation; the remaining 59 tracks! originated in the two-prong stars. In all but three 
of the cases the second prong ended in our part of the stack. Assuming that these 
two prong stars were the products of decay of a neutral particle into a proton and 
a 77 meson, we determined the energy of the decay of each of the 56 such particles 
(see Table I). In 54 cases the values of energy of decay are concentrated in a narrow 


region. 
II. Calibration 


The two quantities to be determined here are the shrinkage factor and the stop- 
ping power of our emulsion. The ranges of 22 protons from the decays at rest were 
measured all in one day on two-Koristka MS2 microscopes fitted with Zeiss Optimeters. 
In order to permit the emulsion to attain a state of approximate equilibrium, an 
effort was made to maintain constant conditions of humidity and temperature in the 
measuring room for about 40 hours preceding the measurement and also during 


the time of the measurement (about 15 hours)?. These conditions were taken as. 


standard conditions to which all measurements of this work were reduced (see$ ZIT). 

In order to avoid coarse errors the results of the proton range measurements 
made on the "standard" day were compared with measurements made by another 
observer and at other periods. Applying to the set of 22 measured proton ranges 
the method developed by Fry and White [2], we obtained for the average value of the 
shrinkage factor: 


sı = 2.20 + 0.033 
and for the mean value of the ranges of the secondary protons 
R, = (1676.2 + 7.3) u. 


The shrinkage factor was also determined by another independent method based 
on the analysis of the set of values Q(s;) calculated by taking the preliminary value 
sı as a first approximation (see § JV, page 281). This procedure gave 


s, = 2.23 + 0.045. 


1 This is less then the number given in an earlier report on this work [1]. In two cases previously 
reported as background, after reexamination, a third track was found. 
? See footnote page 281. 
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The weighted mean of these two independent estimates of the shrinkage factor s, 
and sg: 


2.21 + 0.027 


was used in further calculations. 
As communicated to us by Dr. E. J. Lofgren the density of our emulsion was 
contained between 3.809 and 3.853 g/cm?. Using this information (see [3]) and com- 
bining it with the ratio of the mean range of protons from 27 decays available from 
other work? to our value of R, given above, we obtained the stopping power of our 
emulsion: 
Rs 
R 


This value was used in further calculations. 


— 1.002 + 0.003. 


III. Measurements 


All measurements for each decay were made under high magnification inde- 
pendently by two observers on different days, the mean of the two measurements 
being taken for the calculations. In the range measurements of the secondary particles 
the track was divided into segments which could be approximated by straight lines. 

The measurements of the horizontal projections of the tracks were generally 
made on Zeiss Lumipan miroscopes adapted for emulsion work. The movement 
of the microscope stage could be read on a micrometer screw. In the event that the last 
grain of a track was not in focus with the upper or lower surface of the emulsion we 
made a correction Al estimated from the slope of the track and the distance Ah of 
the last grain from the surface. 

All vertical projections were measured on a Zeiss Optimeter coupled to a Ko- 
ristka MS2 microscope. The measurements were corrected by the formula 


hy = hy ~, - (1) 
x 
where a, denotes the “standard depth", i. e. a previously selected difference in depth 
(~ 200 u of processed emulsion); a, is the same depth recorded the day h, was meas- 
ured. 

The projected angles between the decay prongs were measured on a goniometer 
attached to the eye-piece of the microscope. 

The dip angles of the tracks were measured on a Koristka MS2 microscope. 
In each case the mean dip was determined from a series of measurements of the 


arc tan D ; 


depth h on a base /: 


3 See Appendix. f 
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where the value of / ranged from several to several hundred microns, depending on 
the character of the track. It was found that the variation in thickness of different 
layers of the emulsion as a result of variation of humidity and temperature were 
not strictly proportional. Therefore in the measurement of the dip angles the standard- 
ization (1) was made by means of "partial standards." Three such standards, each 
of a thickness of about 100 u of processed emulsion, were chosen in the vicinity of 
the main “standard.” These standards were used in the dip measurements. 

In order to avoid systematic errors in the ranges owing to the linearization of 
tracks during the measurements, a constant length-correction was added to the mea- 
sured ranges to account for the increased scattering of the track near the final part 
of its range. This correction was determined from measurements of several pion 
and proton tracks, the tracks in question being divided into a large number of small 
‘segments. The mean value of the correction for pions obtained in this way was 10 y, 
while. for protons the correction was negligible. 

An investigation was made to determine whether a small thickness of emulsion 
may have been removed from the upper surface during the processing; this would 
result in loss of part of track. To determine the gap between the upper surface of one 
pellicle and the lower surface of the neighbouring one, measurements were made on 
both surfaces of the distance between prongs from stars located near the surface. 
A number of such measurements indicated that the gap due to air or missing emul- 
sion did not exceed 2 u. In the case of an air gap this factor can be neglected, but in the 
case of missing emulsion it would affect the mean range of steep tracks. À comparison 
of ranges of flat and steep protons from the decays of £7 did not indicate the existence 
of such a difference. 


IV. Evaluation of error. 


On the basis of measurements of the angles and ranges and assuming that the 


secondary particles are protons and pions, we calculated the Q value for each event 
from the formula: 


Q4 = M, — (n, + m,-), 


where 


Mund (Wai Wig)? et Dyce Jes 


W, and W,- are the total energies and p, and p,- are the momenta of ie proton 
atte pion respectively. In the evaluation of the random errors for the individual Q 
values the folloving sources of errors for were taken into account: 

a) Angular measurements. i 

The accuracy of the goniometer for measuring the projected angles was taken 
as 0.5°. In each case comparison was made between this value and the error which 
could be determined from the difference between independent measurements by 


s À 
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two observers‘. The larger of the two values was taken as the error of the projected 
angle. This procedure takes into account both the instrument error and the subjective 
error of the observer. In 26 of the 54 cases the error was greater than 0.5°. A similar 
procedure was used in evaluating the error in the dip angles. The random errors 
of space angles were in 8 cases > 2°, in 3 cases > 3°. 

b) Range measurements and straggling. 

To determine the error in range it was neccesary to determine the error in the 
projected vertical and horizontal segment of the track in each plate. The projected 
length errors were calculated from the differences in the results of two observers, 
separately for mesons and protons; the tracks were divided into classes of large and 
small dip. 

The error in the measurement of the vertical projection was evaluated in a similar 
way (about 5 u per processed plate)’. 

For each track the range straggling error (from [4]) was combined with the meas- 
urement errors. 

The error in each Q value resulting from the above errors in ranges and angles 
will henceforth be referred to as the random error (o, (Q]). 

In the evaluation of the systematic error for individual Q values the following 
sources of error were considered: 

a) Shrinkage factor. 

The shrinkage factor first used for the calculations was s, = 2.20+0.033. Since 
the separation of the A° decays from the background leaves little room for doubt 
(see 8 V) it was possible to determine the shrinkage factor independently under the 
assumption that all the observed Qs refer to a single value. To do this we plotted 


9 
the individual values of Q (sı) as a function of ari Hn is the partial derivative 


taken at the point s, = 2.20 | . In the first-order approximation this dependence 


should be a linear one. From the slope of the straight line drawn through these points 
(weighted with respect to the random errors) we found s, = 2.23+0.045. Finally, 
Lo LC RE ue 08. 7 EM onn Lo RA Eurum o ee S 

4 We calculated the mean error of two measurements by employing Student's distribution 

(3) =18 |x, — #l 
‘ 2 

5 [t turned out that the difference between measurements of thickness of the same emulsion at 
the same place in the plate made under different laboratory conditions of temperature and humidity 
exceeded the expected error of measurement. In seeking the cause for this, several series of measurements 
were made of the emulsion thickness. A number of these measurements were made under different 
laboratory conditions and were standardized according to (1). From the analysis of the measurements 
one may conclude that the “standardizing“ procedure does not guarantee reproducibility of the measu- 


ly varied conditions of humidity and temperature. This introduces 


rements performed under great 1 : 
a source of error considerably larger than the instrument error, the latter in our case being about 1 y. 
ý 


‘ 
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for each case Q(s) and the error ø, {Q} associated with the uncertainty of the shrinkage 
factor was calculated by using the weighted value of the shrinkage factor s = 2.21+ 
+0.027. 


b) Stopping power. 
For all the ranges reduced to the standard emulsion the errors resulting from 


the uncertainty of the stopping power (see $ IT) and the corresponding error for the 
individual Q values (cog {Q}) were calculated. 
c) Range-energy relation. 

Assuming an uncertainty of 0.5 percent® in the range-energy relation [6], we 
calculated the corresponding errors for individual Q values (og (Qj). 
d) Mass values. 

The mass values used in the calculations were [7] 


my = (938.213+0.01) Mev, 


m, = (139.63 +0.06) Mev. 
The coresponding errors in the Q values could be neglected ; however in the calculation 
of the A° mass value this source of error was taken into account. 


V. Results and discussion. 


‘The Q values and their errors for the 56 events mentioned in $/ are listed in 
Table I. The distribution of these Q values is shown in Fig. 1. The Q values of 54 


10 20 30 40 Q in MeV 
Fig. 1 

of our two-prong stars are concentrated in the region 36.5 — 39.5 MeV. In one of 
these events inelastic scattering of the pion with a change of grain density was observed. 
Calculations based on the angle between the prongs, the range of the proton and 
the ionization of the pion gave Q = (39-++7) MeV. The remaining two cases quite 
distinctly lie outside the region 36.5 — 39.5 Mev. The event with Q = 41 MeV occurred 
2 mm from the edge of the emulsion in a strongly distorted region, where meas- 
urements of the angles and ranges are not reliable. The point of decay of anotherone 
(Q — 20.5 Mev) is situated on the very surface of the emulsion or perharps between 
the plates, and consequently could not be carefully examined for the presence of other 
prongs. à 


* Private communication from W. H. Barkas; see also [5] 
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Table 1 
No. Q4 0,{0} 0,{Q0} or{Q} {0} 
1 36.83 0.66 0.27 0.06 0.11 
2 36.85 0.52 0.11 0.07 0.11 
3 37.01 0.44 0.21 0.06 0.10 
4 36.75 0.57 0.20 0.07 0.11 
5 38.36 0.81 0.11 0.07 0.12 
7 38.14 0.62 0.13 0.07 0.11 
8 41.2 
9 37.85 1.24 E02] 0.07 0.12 
10 37.05 1.00 0.24 0.06 0.11 
1 37.33 0.51 0.03 0.07 0.11 
12 36.76 0.58 —0.03 0.06 0.11 
13 37.53 0.46 0.16 0.06 0.11 
14 38.83 0.99 0.30 0.07 0.11 
15 30.22 0.64 0.04 0.07 0.12 
16 37.63 0.51 —0.02 0.07 0.11 
17 3947 
18 37.55 0.67 0.04 0.07 0.11 
19 38.42 0.76 0.25 0.07 0.12 
20 37.99 0.82 0.16 0.07 0.11 
21 38.74 0.56 0.12 0.07 0.11 
22 36.94 0.59 0.02 0.07 0.11 
23 38.61 1.02 0.22 0.07 0.11 
24 20.5 
25 37.45 0.39 0.23 0.07 0.11 
26 37.91 0.71 0.06 0.07 0.12 
27 38.10 0.78 0.02 0.07 0.12 
28 36.63 0.86 —0.08 0.07 0.11 
29 37.87 0.65 0.14 0.07 0.12 
30 37.41 0.62 0.07 0.07 0.11 
31 38.88 0.70 0.07 0.07 0.12 
32 38.20 0.71 0.04 0.07 0.12 
33 37.10 0.50 0.16 0.07 0.11 
34 37.61 0.51 0.08 0.07 0.11 
35 37.08 0.59 0.03 0.07 0.11 
36 37.88 0.47 0.01 0.07 0.11 
37 37.98 0.62 0.07 0.07 0.11 
38 37.60 0.70 0.05 0.07 0.11 
39 37.57 0.52 0.11 0:07 0.11 
40 36.95 0.66 0.28 0.07 0.11 
43 38.36 0.60 0.05 0.07 0.12 
44 37.52 2.14 0.40 0.07 0.11 
45 37.56 0.57 0.06 0.07 0.11 
46 38.80 1.13 0.10 0.07 0.12 
47 38.18 0.61 0.01 0.07 0.11 
48 37.96 0.93 0.23 0.07 0.12 
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Table 1 continued 


7 EL = T 
No. Q4 {0,0} {6,0} ien 9j {080} 

= 1 [ 
49 37.49 0.69 0.08 0.06 0.11 
50 37.88 0.53 0.10 0.07 0.11 
51 37.68 0.56 0.02 0.07 0.11 
52 36.76 0.52 0.11 0.07 0.11 
54 37.12 0.45 0.15 0.07 0.11 
55 38.65 0.54. 0.05 0.07 0.11 
56 37.96 0.53 0.03 0.07 0.11 
1 31.20 0.43 —0.02 0.07 0.11 
58 37.02 0.68 0.13 0.07 0.11 
60 36.50 0.59 0.04. 0.07 0.11 
61 31.88 0.58 —0.02 0.06 0.11 


We therefore conclude that no more than 3 events of our sample are to be regarded 
as background, while the remining 53 correspond to the decay of the AP hyperon. 
The mean of those 53 individual Q values, weighted according to their random errors 
is: 


Q4 = (37.58 + 0.18) MeV. 


The external and internal random errors of the mean are in good agreement: 
0.084 MeV and 0.082 MeV respectively. The final error of the mean includes all the 


random and systematic errors described above namely: 
c, {Q} = 0.08 MeV, 
v, {Q} = 0.09 MeV, 
or{Q} = 0.07 MeV, 
op{Q} = 0.11 MeV. 


À constant area histogram based on the Q values and random errors is shown in 
Fig. 2. 
For the mass of the 4° hyperon we obtain: 
M, = (1115.42 + 0.19) MeV, 


where in the final error the uncertainty in the pion and proton masses are also included 
(c (M4) = 0.06 MeV). Two possible sources of error were not discussed: the influence 
of distortion and the possibility of some undetected inelastic scattering of the secon- 
dary pions. The agreement between the external nad internal errors indicates that 
influence of distortion is small. As regards inelastic scattering of the pions, we note 
that even if one takes the geometric cross section for such interactions (in reality, it is 
probably less [8]), we would expect —2 cases of such scattering in our sample. 
One case of inelastic scattering was found (see 8 V). If we discard even the five lowest 
Q values the mean value of Q} changes by less than 0.1 MeV. 
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area corresponding 
to one event 


E] 


33 36 37 38 39 40 Q in MeV 
Fig. 2 
The evaluation of the 4° mass presented by Alvarez at the Kiev Conference 


on HighEnergy Physics 1959 and based on recent experimental results (which includes 


our value) yields: 
Q, = (37.61 + 0.12) MeV. 


Appendix 


The table below lists data obtained by several laboratories, on the proton range 
from X*— p* 4- 7° decays at rest. All the values are reduced to ranges in standard 
emulsion; the present Warsaw value was based on the density of the emulsion as 
communicated by Dr. E. J. Lofgren. The table gives the total random and systematic 


error. 
Table 2 


Number 


Authors Reference 


of events 


Glasser et al. . 1658.0 8.4 [9] 
Gilbert et al. [10] 
Swami [11] 
Barkas et al. [5] 
Warsaw Group [12] 


Warsaw ‚Group present work 
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The weighted-mean value of the proton ranges according to the data in the table 
above is: 
R, = (1678.6 + 3.2) u. 
Taking into consideration the proton and a? meson masses 
, - 938.213 MeV, 
m,» = 135.04 MeV, 


m 


and the sources of systematic errors: 
A) o{m,} = 0.01 MeV, 
B) c(m,,) = 0.16 MeV, 
C) 0.5 % uncertainty in range-energy relation, 
we obtain the mean value of the mass of the ZT hyperon 
My, = (1189.43 + 0.31) MeV. 
It is worth mentioning that the systematic errors due primarily to (C) exceed 
the statistical error connected with the number of measured decays. 
Om {M,} = 0.01 MeV (A) 
Ons (Mj = 0.09 MeV (B) 
Gp (Mj, =0.27 MeV (C) 


Thus at present the lower limit of the attainable accuracy in the estimation of the 
2* hyperon mass by means of the emulsion technique is 0.29 MeV. 
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ON THE STATISTICS OF SPIN WAVES BY THE 
BETHE METHOD 
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The exact solution of Slater-Bloch's secular equations by the Bethe method with- 
out any additional conditions is considered. The bearing of the results upon the ques- 
tion of statistics of spin waves is discussed. It is shown, that there exists an ambiguity 
which cannot be overcome without additional assumption. 


I. Introduction 


The starting point in the Bloch's spin waves theory is the Heisenberg model of 
a ferromagnetic crystal where the interactions between electrons are considered 
as perturbation. Using the perturbation theory with such a model, one gets in first 
approximation the system of Slater-Bloch's secular equations 


J > {Amy...m, = Qm... } emm E m, m, - (1) 
m.m; : 


The summation is extended over all pairs of contiguous electrons having opposite 
spins, for the indices m,, m; the spins in the pair are interchanged, J is the exchange 
integral between nearest neighbour sites, am,...m, stands for the amplitude of the state 
in which two reversed spins occupy the m,...m, sites and E, is the energy of r reversed 
spin states. | 

The discussion of this system is simplest for a linear chain. For this case Bethe 
(1931) has proposed a treatment in order to obtain a rigorous solution. Basing on this 
treatment Frank (1956) has discussed the symmetry properties of the a amplitudes 
for interchange of k and has arrived the conclusion that for the spin waves the Fermi- 
-Dirac statistics is more apposite than that of Bose-Einstein. Svirskii-Ishmukhametov 
(1957) and Vonsovsky-Svirskii (1959) call in. question this conclusion. They point 
out, the Frank results rest on an artificial introduction of additional equations, which 
requires an extension of the definition of the a amplitudes on the case m, = mj. 

In the present paper this problem will be examined once more without making 


use of any subsidiary equations. 
(289) 
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II. Solution of the secular equations for a linear chain 


The method used in this paper is the same as that of Bethe, i. e. one introduces 
some additional coefficients choosing them in such a way, that the secular equa- 
tions should be rigorously satisfied. The easiest to investigate is the case when there are 
only two spins reversed from the ground state (r — 2). Taking into account the distri- 
bution of the reversed spins in a linear chain, one can divide the system (1) into 
two groups. In the case when the two reversed spins are not located at the nearest 
neighbour sites the secular equations take the form: 


Jan = By tt > Lis ir ann-1 — Anat = Es On (2a) 


On the other hand, when the two reversed spins are nearest neighbours, we get 


J(2a,, 41 mU Amm+2 ve Ct me) = E, Amm+1 : (2b) 
As is well known, the equations (2a) can be satisfied by the substitution 


amni (is ko) = € PR lita en en (3) 


which gives 
E, = 2J {2 — cos k,d — cos kad} (4) 


where d is the lattice spacing, and the k}, kg one sets from the Born-Karman condi- 
tion. The coefficients c, and c, can be freely chosen. One can, however, choose them 
in such a way that (3) will be also a solution of equation (2b). For this purpose we 
substitute (3) and (4) into (25). After a simple calculation, we arrive at the condition 
which the coefficients c must satisfy for (3) to be a solution of the whole system 
(2a) and (2b). This condition has the form: 


. ki— ka 
s (ki E kg) d Ti M TTC 


CR AN ET E 
QR) E 
5 à 


We make now the substitution 


cı (ky, ka) = c3 (ky, ko) 


whence we obtain 
c (k, k) = — ca (k, k) for kAO 
“c; (kys ka) = + ca (ka, ka). 
D this into equation (3) we get (dropping the subscripts) 
amn (ku ki) = c (ky ha) ert L c (kp ko) th Et ka (Ga) 
amn (k, k) =0 k Æ 0. (6b) 
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However, for k, — k, = 0 we obtain from (5) 
AB ky, ko) 
li li Cy ( 1° = 
i0 0 Ca ( ka) 


Taking this into account in (3) one gets 


amn (0,0) Æ 0. (6c) 


These results can be easily generalized to the case r > 2. When no reversed 
spins are nearest neighbours, the secular equations have the form: 


r 
T» (2am,..m,..m, — Qm;..mg—...m, — Am...my+1..m,} = E, Om,...m,...m, + (7a) 
g=1 


On the other hand, when two of the reversed spins are nearest neighbours 
J (2am,... mp, mp1. my. m, — Qm, ...mg 7 1, me 1. mg m, — Qr, .. mp mp 2. mo. m, zn 


r—2 
d 22, dE rom — m,... mp, mg 1. mo— 1. m, — Qi... mp mp1. mg lm, = 
&rs> 


= E, Cm,...mp; mgl. mg m, * (7b) 


(We do not get anything new taking into account more than two neighbouring reversed 
spins). One can satisfy the equations (7a) by the substitutution 


r! "n be 
Am,...m, (ky 500 k,) = Di Cp (ky ee k,) zudem kpj j (8) 
P . 
with 
E, = 2J >) (1 — cos k, d). 
s=1 


In the same manner as in the case r = 2 we look for the condition which must be 
satisfied by the cp if (8) is to be a solution of (7a) and (7b). Substituting (8) and (9) 
into (7b) we get the expression 


rt T x. 

; : ; = 2i k -k 
> cr dE Primi petra + eid paay- epp + ed pan + kp + 
P 


gcfftl 
r TT È 
ipn kp) — D elt pa) — Dy epu ^M) —0 (10) 
ghi z 9 


| ! 
ri i : : 
From all permutations P one can form y pairs of permutations P, P, which satisfy 


the relations 
PO) PTD: 
P,(f +1) =P, (f) | 
P, (g) = P,(g) for all gods (11) 
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This means, that we can group all the elements in > into pairs, which differ only 
P 


in the transpositions of k in the points m, and m, +1. Performing such a grouping 
in (10) we obtain for all pairs 


cp, | cosk— ei? 12 
CP» | COS kext 25, 
where 
p > put + kepu +n) d à — Erein — kpug) d 
cT I D 


Following Bethe and Frank we make the substitution 


Cp = III (+ Y? clkp,, ky) for all P (13) 
i<j 
where the + sign corresponds to the symmetrical and the — sign to the antisymmet- 
rical a amplitudes and p — 0,1 depending on whether the permutation P is even 
or odd. Taking the above considerations into account, we get 


C(kpur, keugty) _ , cosk— e? 


c(kpuct+1)> kp) | ^ cosk— e-* (14) 


Since always one can pick out P, and f in such a way that P, f and P,(f4-1) would 
be any two given numbers from the sequence 1, 2...r, the relation holds for all pairs 


of wave numbers A,...k,. These 3 relations are the conditions under which (8) 


is a solution of the system (7a) and (7b). 


III. Conclusions 


The relations (5) and (14) are deduced directly from the secular equations without 
introducing any additional equations and without extending the definition of the a 
amplitudes. These relations are equivalent to those of Bethe and Frank!. Thus it can 
be seen that the objection of Svirskii-Ishmukhametov and Vonsowsky-Svirskii accor- 
ding to whom Frank results stem from the baseless additional equations, is not valid. 
However, the correct results of Frank do not yet prove the reality of his conclusions. 
The above results are, however, not single valued. In (6a) and in (8) taken together 
with (13) there are namely allowed two signs and the calculations do not answer 
the question which of them is the right one. Since the choice of the sign determines 
RP TN Pe n rockt vor RN 

1 These relations can be obtained by the substitution: 


for Bethe CU E) = eiP(kı,k,)/2 
for Frank c(k,, k,) = e!Mlkı,ka) 
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the symmetry of the a amplitudes for interchage of k, it follows, that the symmetric 
as well as the antisymmetric amplitudes fulfil the system of secular equations. Which 
of them, however, correctly describes the spin waves, cannot be decided without 
any additional assumption. 

This reminds to some degree the situation which occurs in elementary quantum 
mechanics dealing with system containing many particles of the same type. As is 
well known, in this case both the symmetrical and antisymmetrical functions are 
physical admissible solutions of the Schródinger equation and there is no a priori 
way of deciding which one of the two solutions is preferable. This ambiguity is sup- 
pressed by the exclusion principle. 

The application of this principle to the spin waves requires the knowledge which 
is the value of the spin of the quasiparticles corresponding to the spin waves. However, 
this question has not yet been exactly discussed. Generally it is assumed. that this 
value is 1; this would lead to the symmetric function. 

Frank picks out the minus sign in (6a) explaining that it is distinguished by the 
equation (6b) and that for this sign one gets the right number of states in the real k 
space. This argumentation, however, is not convincing because, the equation (6c) 
distinguishes the plus sign. On the other hand, the equation (6b) taken alone does 
not yet mean that the a amplitudes are antisymmetrical. Some light would be thrown 
upon the question by the generalization of the method of this paper for the case 
s >% and consideration of three dimensional lattices. This, however, cannot be 
performed in a simple manner. With more powerful methods this was done by Dyson 
(1956) who came to the conclusion that the quasiparticles of spin waves can be con- 
sidered as a Bose gas, subject to no exclusion principle. The investigation on the neutron 
scattering, performed by Riste, Blinowski and Janik (1959) seems to support this 
conclusion. 

I am indebted to Prof. S. Szezeniowski for his continuous interest and encour- 
agement during this work as well as for reading and discussing the text. I wish 
alsa to thank Dr. H. Cofta for discussion on some points. 
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EFFECTIVE MASS METHOD IN THE CASE OF 
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The case of a non-parabolic energy band in a semiconductor has been analyzed 
and the necessary modification of the usual effective mass method has been found. It 
turns out that to describe the motion of single electrons, one must introduce instead 
of the coefficient mef several coefficients, called the differential effective masses. It has 
been shown that in the above case-the phenomena of electrical transport can be described 
by the same expressions as in the parabolic case with the only difference that m,; has 
to be replaced by a characteristic quantity for the given phenomenon, the so-called inte- 
gral effective mass, being a mean of the differential effective masses. Making use of E. O. 
Kane's dispersion formula some integral effective masses for In-Sb — n have been cal- 
culated, as well as some scattering coefficients appearing in the Hall and Seebeck effect. 


The theoretical description of transport phenomena in semi-conductors is usually 
based in the so-called effective mass method wherein the carriers in the semi-con- 
ductor are assumed to possess the properties of free electrons the mass of which, 
however, has undergone modification. It is this modified value of the mass (effective 
mass, m) that accounts for interactions between the lattice and the carrier. . 

Thus, the effective mass method represents a problem of mass renormalisation. 
The method itself was developed on two tacit assumptions, namely: 


4 


im 


The assumption of spherical structure of the conductivity band and of the fun- 


damental band, as expressed by the condition € (E) = e(k), with e denoting the energy 
of the electron, and i. its wave vector. Until recently (up to the date of discovery of 
the cyclotron effect in Ge) it was considered to be justified in the case of all crystals 
of the regular system. As a great many semi-conductors and, first among these, those 
that can be obtained to a very high degree of purity, crystalize in the regular system, 
the assumption was thought to be fully justified. 

(295) 
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Il. 


Moreover, the electrons (holes) were assumed to occupe only levels sufficiently 
close to the bottom of the conductivity band (top of the fundamental band) for € 
to be given by the first term in the Taylor expansion of e(k) in powers of k. Clearly, 
at the extremities of the band, the Taylor expansion of the energy begins from 
the quadratic term; thus by assumption (I): 


E = 5 pa ne k? + higher order terms (1) 
or, denoting 
i(s 2r 
h? 9r? on Mef 
and with hk = p. wherein h = h = Plancks constant 
Sle 
a= m (2) 


Assumption (II) was based on the fact that the free carriers usually occupied a very 
small part of the band in semi-conductors (i. e. substances wherein the permitted 
bands are separated by energies of the order of 1 eV). Even at medium concentrations 
of the order of 1016 cnr? only one millionth percent of the free states is occupied. 
Hence, in a great number of cases, this may be a region wherein the quadratic ap- 
proximation is really sufficiently accurate. 

With the foregoing assumptions, the behaviour of the electrons in the crystal 
when subjected to an external force F is described, according to Bloch's general 
theorem (1929), by the classical law of motion: 

A E LE (3) 
(with ? denoting the mean velocity of the electron in the crystal), the only difference 
consisting in the fact that instead of the rest mass of the electron, eq. (3) contains 
the effective mass of eq. (2). 

On the other hand, eq. (2) yields a simple expression for the density of electron 
states (the number of electron states per unit energy interval): 


D (e) = F Qm, VE (4) 


The simple effective mass method based on assumptions (I) and (II) will be 
here termed the "classical" method, with respect to its close bearing on the classi- 


cal theory of electron motion and the part played by it in the oleis. of 
semi-conductor physics. 


- 
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It will be subsequently proved that the inverse is also exact, namely: the assump- 
tion of constant carrier mass leads to the conclusion that assupmtions (I) and (II) 
must be fulfilled simultaneously. Thus, “classical” effective carrier mass implies 
a highly specialized (simple) energy structure model of the semi-conductor. Inverse- 
ly, a substance for which Assumption (I) or (II) is not fulfilled is not accounted 
for by the “classical” concept of effective mass. 


Renouncement of Assumption of Sphericality 


Investigation of the cyclotron effect in Ge and Si carried out in 1954-1955 
(A. F. Kipp, C. Kittel 1955) clearly proved that Assumption (I) (that of sphericality) 
was not justified in the case of these crystals. Experiments showed the conductivity 
band to be given by a dispersion formula, elh), corresponding to surfaces of equal 
energy having the shape of greatly elongated rotational ellipsoids whose principal 
axis had the direction of [111] (Ge) or [100] (Si). Here, minimum energy does not 
coincide with k — 0, but with several crystallographically equivalent points k0 
wherein i = 1, 1, ... N, and N denotes the number of minima of e(k). 

Hence, this represents a model that is not only non-spherical, but is at the some 
time one possessing several minima points (multi-valey model) In the neigh- 


bourhood of each minimum in point A, the dispersion formula is of the shape 


h? ji — ki)? , (i5 — ki)? + (ks — nA 
Mem n EE ns 


m m 


e$) = (5) 
In eq. (5), the principal axes of the i-th ellipsoid have been chosen as coordinate axes. 
Clearly, the coefficients m’ and m" have the dimension of mass and replace the effec- 
tive mass of earlier formulas. (see eq. (2)); the quantities m’ and m" are usually termed 
the longitudinal and transversal effective mass, respectively. 

Eq. (5) can be rewritten to become formally identical with eq. (2) viz.: 


wherein Pi =h (k = $’) (2°) 


The constant quantity m,, of eq. (2) is now replaced by one that is, generally, 
variable: m* — m*(k). Considering eq. (5) to represent the specialized dispersion 
formula, it is readily seen that m* depends solely on the angle subtended by the 
vector k and the rotation axes of the ellipsoid; thus, on changing from the rectangular 
coordinates k,, ka, k to the spherical coordinates k, 9, p, the quantity m* is a function 
of the polar angle 9, thus: m* = m* (0). The dispertion formula is clearly seen to 
conserve the quadratic form analytically, i. e. notwithstanding the fact that the first 
condition of the former model has been removed, the second condition is still ful- 
filled. The clasical equation (3) will be replaced by the law of motion for a Bloch 
wave packet in the presence of.external forces, in tensor form (F. Bloch 1929). 
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Thus, in place of eq. (3), we have now 


cr (3°) 


wherein the tensor E (in the system of principal axes of the ellipsoid (5), as above) 


is given simply in the form 


"ou 
m 
1 1 (6) 
x c Saa 
emen 
m 


Thus, eqs. (29) and (39) replace eqs. (2) and (3) of the common effective mass theory 
in the case under consideration. Eqs. (29) and (39) still retain the form of dynamic 
laws and can be used for investigating electron motion in a crystal by classical me- 
thods (common diffrential equations). The present case, however, introduces a more 
profound modification into the laws of motion of the electron than that implied by 
the effective mass theory in its "classical" form. 

In the present case, both fundamental laws of motion (29) and (39) no longer 
contain one and the same effecive mass, but two distinct “mass coefficients" namely: 
the energy formula contains the quantity m* = m* (m', m"; 9), whereas the force- 
acceleration relation contains no less than 9 tensor elements which are yet other 
functions of the same coefficients m’, m" and 9. Characteristically for the case under 


: : 1 , 
consideration, both m* and the tensor — are independent of the absolute value 


M 


of the wave vector k but depend solely on the angle subtended by k and the rotation 
axis of the ellipsoid. The foregoing fact, which is a consequence of the quadratic 
form of eq. (5) implies quite generally, in the expressions accounting for various 


: 1 
electronic processes, that the structural factors related to the quantities m* and 4; 


are independent of the degree of ivestment of the band (and, hence, are concentra- 
tion and temperature independent). 

The foregoing relatively simple example of failure of the former concept of effec- 
tive mass has been dealt with at length in order to show, that — if the formalism > 
is appropriately modified — the motion of the carriers, and, thus, a specialized trans- 
port phenomenon can still be described by a method that is, in principal, classical. 
On the other hand, the computations required in the case under consideration as 
given by eq. (5) for the electronic effects hitherto accounted for by the “classical” 
method of constant effective mass are, as yet, a highly involved problem and one 
whose greater part still awaits solution. It is not the present author’s intention to 
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give a detailed account of the computations, wich will be reproduced here in their 
general outline only. Thus, eqs. (29) and (39) should be completed by an expression 
accounting for the density of states, formally idencitcal with eq. (4), but such that 


mis therein replaced by the quantity m; being yet another, though well-defined 
function of the coefficients by eq. (5), namely 


4 
m, Nim tm * (4°) 


To replace (2), (3) and (4), the quantities (2°), (3°) and (4°), respectively, should be 
substituted in the expressions for the various electric effects in semi-conductors. 


] 
Hence, expressions with m*, M and m,, will appear in the respective sums and 


intergrals resulting from integration of the momenta or energies over all the electrons. 
This means that the formulas for the transport phenomena retain, in the case of Ge, 
their previous formal character, with various new structural expressions (related to 


the form of e(k)) to replace in the various formulas. The generally different structural 
expressions will be jointly termed mass coefficients or effective masses, each of these 
being additionally denoted according to the kind of phenomenon under consider- 
ation. In this manner, the universal effective mass of the simple model is now replaced 
by a set of different masses. All mass coefficients are obtained from the three funda- 


mental masses m*, Me "as. and can be ultimately and univocally derived from 


the dispersion formula by the procedure already outlined. Here, as exemplified by 
the ellipsoidal modification of the band, a consistent procedure has been put to effect 
making it possible to correct the theory. This requires, in general, that all three fun- 
damental formulas defining the effective mass be modified. 

In a paper on germanium (Z. Kopeé 1956), the formula accounting for the ther- 
moelectric force and mobility has been corrected along lines corresponding to the 
modified theory already mentioned. The corrected expressions have served to remove 
certain experimental “anomalies” existing in this field of research and have made 
it possible to relate quantitatively the various electronic phenomena in germanium, 
thus assembling them into one consistent group. In particular, the (by “classical” 
standards) anomalously high thermoelectric force in Ge — n as compared with Ge — p 
found its explanation; similarly, the source of the discrepanéy between the effective 
mass values from electron mobility, on the one hand, and thermoelectric force, on 
the other, was clarified (the respective values had differed by 150—200 %). Within 
the framework of the corrected theory, the foregoing “anomaly” served to corroborate 
the discovery of band anisotropy. 

In the analysis of carrier effective mass in Ge as proposed at the time, the author 
based on eq. (5). It has been already mentioned that this is equivalent to retaining ` 
Assumption (II) of the previous concept of mass. 

Sufficiently far from the bottom of the band, however, there should be divergence 
from the quadratic law for £ = f(k), leading to temperature and concentration depend- 


300 Z. Kopeé 


ence of the effecitve mass. The experimental results obtained with Ge — n samples 
of a very high degree of purity are in agreement with computations based on the 
quadratic approximation. On the other hand, in samples of medium and high carrier 
concentrations, numerous measurements lead to disagreement with theoretical results 
based on this approximation (Z. Kopeé 1956), which points to the probability of a di- 
vergence of the kind discussed. However, at the present stage of our knowledge 
of band structure in Ge — n, no quantitative expression of this divergence could 
be given. 

In InSb—n, the state of affairs is one of still much more drastic desagreement. 
In this case the effective mass of electron state density is by over one order lower 
than in Ge. Hence, levels that are far remote from the extremity of the band are 
occupied already at relatively low concentrations. Thus, we are now entitled to expect 
especially grave results of the divergence from the quadratic dispersion formula. 
Indeed, such divergences were shown to exist from our experiments on the thermo- 
electric force in InSb and from the optical experiments by Fan and Spitzer on infrared 
absoption in strongly doped InSb samples (H. J. Fan, W. G. Spitzer 1955). 


Renouncement of Assumption of Quadraticity 


To generalize the effective mass method to include non-quadratic structures 
(eÆak?), a quite general problem will now be discussed. No assumptions whatsoever 
will be made as to the form of the dependence of the energy of the electron on the 
set of quantum numbers accounting for its state. For a dependence of quite general form 
it will be the author's aim to derive laws of motion formally analogous to the 
entirely classical formulas, and expressions for the density of states formally identical 
with that of the “classical” theory of effective mass discussed in the Introductory 
chapter. 


In the equations describing the motion of the carrier in external fields and in 
the fields of various perturbations of the ideal crystal lattice (to which energy band 
structure relates) different effective masses will appear which in general, can be 
functions of all the quantum numbers accouting for the state of the electron. Th 
us each electron will in general, be described by several distinet masses: the energy 
mass, state density and dynamic masses. Each of these quantities can be different for 
the various electrons within one and the same semi-conductor. 


By summation (statistics) of elementary processes, it is possible, as previously, 
to achieve description of phenomena involving great numbers of electrons (trans- 
port and other electronic processes). The structural factors appearing in this de- 
scription will already be highly complex. In following we find namely the form of these 
coefficients in nonparabolical case. 


All mass coefficients can be classified in one of the following two groups: that 
of differential and that of integral mass coefficients. | 
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Differential and Integral Effective Masses 


Quantities accounting for a separate electron in the band will be termed diffe- 
rential quantities. Three fundamental effective masses will be distinguished, in accor- 
dance with the following three definitions: | 


P 
TS (29 
l dv; 
et I 
m" 99 
4 
Die) = s (2ma.. ^ Ve (4) 


YE 1 
The quantities m* and m account for the properties of motion of the electrons 


in the semi-conductors in a field of external forces; on the other hand, m4, de- 
scribes the electron state density D(e) whose knowledge is indispensable for summation 
over all occupied electron states in the band in the procedure of computing transport 
processes. 

Eqs. (21), (3!) and (4!) can be considered to arise from (2), (3) and (4) by replacing 
m,, therein by m*, M and m,,, respectively. 

To conserve the previous form of the expressions describing transport phenomena, 
integral effective masses, being certain combinations of the differential masses, will 
be introduced. 

In the general case, any electron transport phenomenon was accounted for in the 
classical effective mass method by relation (a): 

(a) z — z (semi-conductor parameters, external conditions, m,,.); in the present 
case, in accordance with the generalized fundamental equations (21), (3!) and (4!) 
the expression for z will be appropriately modified. This is carried out as follows. 
In a great number of cases, the expression corresponding to the phenomenon under 
consideration will be able to be put in the form 

(a) z — z (semi-conductor parameters, external conditions m,,); identical with 
the “classical” form, but one wherein m,, has replaced the previous universal effecitve 
mass m, . For a given phenomenon, 7j will be expressed by a specific mean value 
of the differential effective masses, to be termed the integral effective mass of the 
electron phenomenon under consideration (index: p). 

In some cases, a structure-independent expression can correspond “classicaly”” 
to an electron process: 

(b) z' = 2’ (semi-conductor parameters, external conditions). 

On the other hand, within the framework of the generalized effective mass method, 
this kind of structural dependence can occur, then assuming the form 


LÀ 
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(b) z =z (semi-conductor parameters, external conditions; 2j). wherein Ap, 
represents a structural coefficient being an integral function of the same differen- 
tial massess. 

The possibility of the foregoing simple generalization of the "classical" effective 
mass method results from the fact that the “classical” quantities (2), (3) and (4) are 
replaced by the three fundamental equations (21), (3!) and (41). 

The set of mass coefficients play a role similar to that of the effective mass of the 
simple theory and lend profounder significance to the intuitions inherent in the 
*classical" method of effective mass. 

We shall now proceed to compute the differential and, subsequently, some of the 


integral mass coefficients for an InSb crystal. 


Differential Effective Masses in InSb—n 


Using the wave vector as quantum index of the state of the electron (perfect 
periodicity), the equations can be written as follows (see Appendix): 


h2k2 


en 9I 
E— 8 (211) 
ETC AL: IH 
M; h? 3k;3k; Se 

1 h? -D (k, ke, ks) i 
= |e E dd 4n 
pe las Dereng ^ * idein 


It should be noted that eqs. (2H) (311) and (4) define the fundamental differential 
effective masses in accordance with the scheme assumed for a perfect crystal (de- 
scription in terms of the wave vector). However, eqs. (21), (37) and (4!) are of a much 
more general form. Eqs. (2%), (3!) and (4%) result from definitions (27), (3!) and (45) 
as applied to the Brillouin band (description in terms of the wave vector h). However, 
various experimentally known facts speaking in favour ofa qualitative analogy between 
phenomena in periodical and amorphous structures (vitrous semi-conductors) seem 
to indicate that the possibilities of a semi-classical description of the dynamic proper- 
ties of the electron in solids are not restricted to conditions of perfect periodicity. 
For non-periodic structures, eqs. (21), (3!) and (4!) can be assumed to represent original 
definitions of the quantities m*, = and mg, whereas the state of the electron 
can be described by an adequate set of quantum numbers in place the wave vector k 
characteristic of the state of the electron in periodic structures. We shall now proceed 
to compute the values of the differential effective masses for electrons in InSb. Here, 
a number of experimental results point to spherical structure of the conductivity 
band, (thus, e. g. experiments on variations in the electric conductivity in a longitu- 
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dinally aplied magnetic field) Hence, the dispersion formula for InSb will be of the 
form 


e(E) = «() (7) 


With respect to the non-quadratic form of e(k), it would be possible to take 
into account the next, higher term in the Taylor expansion, thus assuming 


elk) = ak?-- bk? (8) 


In the next step, with this expression, it would be possible to determine various 
effective mass quantities and to find the values of a and b from comparison with 
experimental data. If necessary, powers higher than 4° in the expansion of e(k) should 
be taken into account. 

Here, however, a different procedure will be adopted. Indeed our starting point 
will be given by Kane's formula (9), (1957), recently derived along theoretical lines, 
which explicitly gives e(k) in InSb. 

Hence, the exprimental confirmation. of the results of the present computations 
can be considered to provide the justfication of our general procedure and, at the 
same time that of Kane's formula. 

Kane proposes the following expression for the energy of an electron in the 
conductivity band: 


—h 1 D N | 
e (k) UU RO a + B h2k2P Aec (9) 


The formula contains, in addition to universal constants, only two structural para- 
meters which should be deterrnined experimentally, namely : 

Aeg the width of the energy gap, and the quantity P, which can be related to 
the cyclotron effect and determined from formula (10): 
with m, denoting the cyclotron mass (see Appendix). 
Zp = Aeg — 


Mo * Me 


— Me 


(10) 


Eq. (9) was derived by highly interesting consideration of a general nature not 
involving knowledge of the form of the InSb lattice potential, but basing on its sym- 
metry properties and characteristics as resulting from well-known and rather simple 
experimental facts (a single resonance maximum in the cyclotron effect, behauiour 
of the electric resistivity in a magnetic field, an extremely narrow forbidden band). 
Thus, Kane could adopt the most appropriate perturbation method and restrict the 
problem to kp and spin-orbital interactions indepedently of k. Kane’s formula ınay 
indeed be considered to account adequately for the structure of the conductivity 
band, at least for electrons of up to 0.3 ev, since, as that author proved, neglected 
| interactions of various other kinds begin to play a part for electrons whose energies 


exceed this value. 
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<] 
The energy mass, m*, as determined from Kane’s formula, is conveniently 
represented in the form of an expansion. It is readily proved, from eq. (9), that B 


emeret eg" eset 


(11) 

The differential state density effective mass m,,, is expressed in m* according to 
d © * "ls i i 

mas, = m*h (m + °) | ineat (12) 


1 
The tensor M is given by the following expression: 
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It is noteworthy that, notwithstanding the fact that the band structure is sphe- 
rical, the effective mass tensor M by no means reduces to a scalar. Hence, the 


change in velocity of the electron in the field of force does not, in general, oecur 


SIR 


'06 


m, 


"05 


04 


/N. 
2 
Mas. 


03 M. 


02 


01 


.2 E à ‘8 10 E/sEg 


Some differential effective masses for InSb. Effective energy mass — m*; state-densities — mq; dyna- 
mical masses as functions of electron anergy (calculated with m, = 0.014) — mi, mg. 


in the direction of the force, but depends also on the direction of cristal-momentum 
of the electron. l 

Scalarity of the effective dynamic mass 1s exclusively a property of a structure 
that is at the same time parabolic and spherical (see Appendix):, it thus represents 
a case that does not occur in normal conditions among the semi-conductors known 


at present. 
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Some Integral Effective Mass Quantities in InSb—n 


State density mass Mas.. 
The quantity 1s defined by the well- known formula accounting for the electron 
concentration in the band: 


2(2mkT)h Ma, 2 
LE h3 ya Fy, (n) (14) 


wherein 
k — Boltzmanns constant 


h — Plancks constant 


ne f- and g — energy of Fermi level 


kT 


~ 


F; (n) = -f i LES (Fermi's integral) 


0 

The expression can also be derived in the Kane’s type band. The quantity Mas. 
is obtained from the differential coefficient m,, by integration according to eq. (15) 
(see Appendix): 


e "Is 
fo ma, hs x 
1 + exp t = PRU (15) 
Fy, (n) 


Exact calculation of M,, from eq. (15) reduced to the numerical computation of 
a sequence of Fermi type integrals and is a highly toilsome affair. Since, however, 
the &(k) dependence in not expressed quite exactly by Kane’s formula (especially 
at energies € > 0.3 eV), it seems justified to introduce an approximation that greatly 
simplified the computations, albeit at the cost of some small increase in error. 

It is readily proved (see Appendix) that, assuming the linear approximation 


Mass == 


sh=s+t (16) 


"ms 


instead of the exact expression for m, ^in eq. (15), a relatively Simple, approximate 
form of M,, is obtained: 

kT Fy,(n) 
Aeg Fy(n) 


The following values have been assumed: s = 0.00165, t = 0.0055. The error of 
this computation does not exced 10 %. Here, the values of m, and Aeg had to be 


LEE Ter (17) 
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taken from other sources. For m, ,a value of 0.014 (in my) was assumed, representing 
the average of m,!) as measured in the cyclotron effect at 4.29 K: m, = 0.013, and 
at 300° K: m, = 0.015. 
dA£g 
dt 

in accordance with the optical data found in the literature. 

Some results of the computations are given in Table 1, wherein the coordinates 
T(°K) and n (concentration of free electrons) are employed, as is usually the case 
in describing the state of the electrons in a sample. 


The value of Aeg (s OK Aeg = 0.24 ev —217: tote was assumed 


Table 1 Mj, (T, n) 


PES] 1015 1016 1017 1018 1019 


.0148 .0151 .0172 .0239 .0394 
.0155 .0159 .0176 .0242 .0410 
.0194 .0198 .0206 .0266 .0451 


From Table 1, the effective state density mass in InSb increases with the temperature 
and the electron concentration. The author feels he should stress that, at 3000 K. 
the transition from electron concentration 101$ cnr? to 101? cm”? involves an increase 
in M,, to more than double. It is this integral effective state density mass that is 
obtained in measurements of the thermoelectromotive force œ (see Appendix). 


Dispersion Coefficients A and r 


In the expression of the thermoelectromotive force 


- (4-5 (18) 
where e denotes the elementary charge, the well-known dispersion coefficient À repre- 
senting the mean kinetic energy of the carriers in kT units appears..The coefficient 
itself has been the object of many an investigation, whence it is known to possess 
a number of interesting properties. 

In the first place, for applicability of Boltzmann statistics and of the parabolic 


band approximation the values of À for most dispersion mechanisms range from 2 


1) If the result of measurements of m, in the cyclotron effect at room temperature is analyzed 
more closely the conclusion is reached that the true value of the effective mass at the bottom of the 
‘band at the temperature considered is somewhat smaller than the value of 0.013 characteristic of low 
temperatures. Since the computations were approximate, the present author refrained from introducing 


the temperature variations in mç Besides, to obtain the latter variations, some kind of extrapolation 


would have to be adopted. 


4 
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to 4. A value of 2 corresponds to dispersion on acoustical phonons, whilst that of 4 
to dispersion on ions. If both mechanisms exist simultaneously, intermediate values 
appear; the dependence of the latter on the rate of participation of either mechanism 
is discussed e. g. in the paper by A. I. Anselm and V. I. Klatchkin (1952). Degeneracy 
has been investigated by a rigorous treatment only for pure dispersion mechanisms 
(J. Tauc 1953, O. Madelung and H. Weiss 1954). In cases when dispersion can be 
accounted for by introducing the notion of relaxation time, the formula for 4 assumes 
the form 


4x mee (19) 


with 


soe ee inia Do , 
ls - $e f» l e D (e) 5— de (20) 


The following relation can be proved to hold for dispersion on acoustical phonons 
and Kane’s structure (see Appendix): 


7 hi Cil m* 


OT ma 


(21) 


being a modification of the respective formula of deformation potential theory (J. 
Bardeen, W. Schockley 1950), where as D(e) is determined from eqs. (4!) and 12); 
finally, 4 is obtained in the form 


URL PE LR." (22) 


wherein 


x= elkT and ET 


RON 


For strong degeneracy, with respect to the deltoidal form of ps A — n. The differ- 
x 


ence (A — n), being a second order quantity, requires a separate calculation. In the 
present case, the result will differ but negligibly from the one obtained for the parabolic 
band. 

The situation is entirely different in the absence of strong degeneracy. For Bolt- 
zmann statistics, we have 


= PRATER AH (221) 
"T xe-* dx 
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+2 

m 

at energy values of e < 0.3 ev, the quantity MS can be approximated as 
S 


follows {at higher energies, the contribution to the integral in negligible): 


m*2 


= ao +a, x + a4 x? (23) 


m, ? 


the quotient of integrals (22°) reduces to the approximate expression 


__2!a+3!a+4læ 
ASISTE. ER 


assuming m, = 0.014 and Aeg = 0.175 ev; A = 1.66 (T = 300? K). 

If dispersion is not phonon acoustical but phonon-optical, the quantity 4 can no 
longer be expressed by eq. (22) (as it is not possible to define the quantity l, and 
the problem becomes highly involved nathematically; Boltzmann's equation then 
becomes Voterra's integral equation, which must be solved by variational methods). 
A case of this kind was discussed by E. H. Sondheimer and O. I. Hogath (1953) and 
applied to Kane's structure by H. Ehrenreich (1957). According to the latter author, 
A = 2.5 at 300° K. 

As to the coefficient r in the formula 


r 


me (24) 
with : 
R denoting Hall's constant 
n — the carrier concentration, 
c — the velocity of light, and 
e — the elementary charge, 


it is independent of the dispersion mechanism at high field intensities, and assumes 
the value r — 1. On the other hand, however, non-parabolicity of the structure, 
by a computation analogous to the one yielding À, leads to an increase in run enl: 
For strong degeneracy r again tends to 1. 


Phonon -Acoustical Mobility Mass 


We shall also deal with the mass coefficient in Shockley's well-known mobility 
formula: 


|. (8n) ^ eh* cu 


Ace 25 
Koh mls (ET)! e? es) 


It is the general opinion that, in pure germanium and in pure silicium, the car- 
riers at room temperature are dispersed on the acoustical phonons exclusively. 
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Consideration of the steps leading to 7 in eq. (29) yields a relation between 


: i ; - a 
m, and the integral state density mass M,, already derived (in the case of Kane's 


band structure — see Appendix): 

my = Ube m, ^ (26) 
The quantity m, in the foregoing formula is given by eq. (27): 

2 2 

mn AO) Are vi) 


Me = © (27) 


"mas exp (x — n) 1 
en MES | 
mas? {1+ exp (x— 7)} 


We, = 


It is seen that m, is closely related to the coefficient x, in eq. (22) for 4. From eq. 
(27), m, was computed by expanding in a series of Fermi functions and taking several 
terms, with an accuracy of 15 — 20 %. Substituting the respective values of m, 
together with those of M,, from Table 1 in eq. (26), the dependence of m, on the 
concentration n and temperature Ÿ was obtained. The results are shown in Table 2. 


Table 2 m u 


1015 1016 1017 102 


.0156 .018 
.0165 .018 
.0193 .025 


The vaıue of m, is seen to exceed considerably both the effective cyclotron mass 
and the effective state density mass. The concentrational variability of the quantity 
considered is a much more pronounced one. It is highly characteristic that, contrary 
to M,,, the quantity m, can be a decreasing function of the temperature (at constant 
concentration). Namely, in sufficiently impure samples, m, decreases with higher 
temperatures. ! 

The foregoing results of the computation of mass coefficients in InSb, and, 
more especially, the Table of integral effective state density mass values (Table 1) 
and mobility mass values, point not only to far reaching divergences in the numerical 
values of these coefficients, but also to a different form of their temperature and 
concentration dependence. 

These are the forms — certainly “anomalous“ ones from the classical point 
of view — that have been observed iu investigations on the thermoelectromotive 
force and carier mobility in InSb. An account of these experiments will be given 
elsewhere. 

The author wishes to express his very great indebtedness to Proffessor A. F. 
Joffe for his valuable discussions. It was Professor Joffe's initiative that the author 
took up the problem of effective mass in InSb. 
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Appendix 
State Density 


In the “classical theory of effective mass, the density of states (i. e. the number 
of states per unit energy e and ünit body angle 2) is given by the relation 


De We 9 
d ys (2m) y2 


d (4a) 


which is independent of the angle (isotropy); hence, it is always possible to employ 
state density per unit energy: 


D (c) = 4nd (e) = P. (2m, Ve (4) 


Quite generally, in the description of the state of an electron by a set of numbers 
Xo Xy X in x-space, a function dọ (x, ...) will also exist, which we shall term the state 
density per x-space element: 


do(x4, Xos x9) dv, = aN 


with dr, = dx, dx, da, and dN = the number of states per element dr,. 

In considering an isotropic problem (the only case reasonably possible in genera- 
lizing over homogeneous and amorphic structures), d,(x) will be a function of the 
energy only: 


da) = de) - (4b 


We shall now write an expression, analogous to (4a), for the density of states in the 


general case: 


er (2m,.,.)'^ Ve (4a) 


do is 


which we shall consider a definition of the generalized differential effective state density 
mass My, 

In concrete cases of computation, we are interested in periodic structures only. 
Isotropic periodic structures (InSb) are also entirely within the scope of the foregoing 
general notation (4). Anisotropic structures, however, require that we start from 
the more fudamental equation (4a). Clearly, the differential state density mass, as 


defined by eq. (4a) can now be a function not only of the energy, but also of the angle: 


"a 
d(e, 0,9) = 2 E (4a) 


The form of m,,, as defined by (4a) is obtained through finding the Jacobian 
of the transition from the coordinates ky» ka kg to £, 9, y. In the former system, the 
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state density in an element dt, = dk, dk, dk, is constant and amounts to Qa 


9 
ann y = 4 
am IN (4c) 
We have 
dk, dk, dk, = Dea o) de dà dg xis Dia p de (4d) 
Substituting (4d) in (4c) we have 
2 1 D (k kaka) 
D (e, 9, 9) = de: dQ 4 
(22)? sind? D (e, 0, y) de dQ = d (e, 9, y) de (4e) 
1. ©. 
2.1 T 
~ (2x)? sind D(s,0,9) 
or 


D (ky ka ka) 
M C De ^" 4 
and, by the difinition of eq. (41) 
D (k; ka ka) p|" 
; ddd qu 
ne z lay: D (e, 9. 9) d 


Kane’s Formula 


In numerical computations, we employed the following form of Kane’s formula: 


ERAS EE 


= Om, 


Ya 
242 UFKA EE I 
+4 | deb + 20242 400 Im me) Ac | (91) 
Eq. (9) is derived from (9) by substituting eq. (10) therein. The quantity m, defined 
in (10) is the cyclotron mass, which can be obtained from measurements at low tem- 
peratures with samples of very high purity. 
In these conditions, it is the electrons at the very bottom of the band that come 


to be investigated, for which k 0. By transition to the limit in eq. (9) for k>0 
it is readily found that 


h?k? 


E > 
2m, 


ko — 


(9a) 


Thus, for such electrons, we have once more derived the quadratic approximation 


Hence, in their case, the “classical“ concept of effective mass holds and in cyclotron 
measurements Mef = Me: l 
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Effective Mass Quantities in InSb 


From the definition (2") and eq. (91), the energy effective mass m* can be deter- 
mined in accordance with eq. (11). 


€ 
The form of eq. (11) which is that of a simple power series in powers of EE 
G 
is well adapted to the problem of determining all the remaining effective mass quan- 
titiens in InSb as the latter will be given by the expansion (11) together with its first 
and second derivatives: 


STER TENE el On m3 5,1 ds ms Capuc ee 
ee DEK 1) n! RER: » (me) Mo‘ Me Aeg 


(lla) 
©?m* 
9 (ej As)? IS 
= i Me n—2 [m S d n—3 
=m y (—1 EA) rn) (n —2) 6.) poen be) 
n=1 
(11b) 


By simple computations, in accordance with eq. (41) and with respect to 


D (ky ka ka) _ 2m*e 9k 
D(e, 9,9)  h? de 


wherein 


1 es Om* 
9k P m* Aeg 9e[Aec 


= T 
de hAec'h | DE 


m* Aeg 


the differential state density mass reduces to 


Om* € dire 
Ye Eee eee i 12 
HUG ie Im" Im + SE PS | ( ) 


Eq. (13) results immediately from Bloch’s theorem on the motion of a wave packet, 


as applied to a band described by a formula of the form e(k) = e(k). It is seen that 
the reciprocal mass tensor can also be computed from eqs. (11), (11a) and (11b). 
By the form of eq. (13), the dynamic mass is seen not to be a scalar, to the sole 


exception of the case 


3e 1 de 
DE kə 


i.e. € = const. k? 
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Hence, the parabolical and spherical energy structure, i. e. the case of the “classical” 
concept of effective mass, is the only one to admit of the scalar form to the effective 


dynamic mass. 


Integral Mass Quantities 
State Density Mass M,,. 


On integrating the product of the state density D(e), as given by eq. (4) and 


of the electron distribution function (Fermi) f(e) = over the 
1 + exp Ca ) 
energy, we obtain the number of electrons in the band, n = fD(e) f(e) de. 
On expressing D(e) by m4, according to eq. (4), we have 


[E m 

2(27 kTyh 2 1+ exp (x — n) 

= Fy, (yn) I — —— (14) 
7L h3 Va la ( 1) Fy, (n) 


Thus, it has been possible to express the concentration by an equation that is formally 
identical with the “classical“ expression (the “classical“ concept of effective mass) 
with the only difference that the respective power of the effective mass is now replaced 
by the expression in square brackets. Denoting the latter, in accordance with the 
general procedure adopted in the present paper, by [ ] = M, ^ the definition of the 
integral state density mass is obtained: 


f mas. xh fda]’» 
Mem [Rey Le 


Dispersion Coefficients 
Here, the non-quadratic form of e(k) should be acounted for in the expression 


for the mean free path J or the relaxation time v. According to W. Shockley (“Elec- 
trons and Holes in Semi-Conductors“, 1950), 


1 
E EA J Wo; 42, (D) 
wherein E 
: Uol? 
Wo = RL ono; (1) 


the probability for transition of an electron of energy €, from state Ó to state j, |Ujo |? 
denoting the square of the transition matrix element (Uo lE = en o; is given by the 


relation: 


V b. bs > 
73 di dps dps = oj dQ; de (P= hk) | (III) 
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in the case under consideration 


dp dps dp = " men dO;de IV 
as 
| pel = v, = Dem (V) 
E 
dp, dp,dps — 5 dQ; de (VI) 


V Vernier EMT. 
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D9 o a a ———————— zz 


*2 *2 

m m dom 
Approximating —— and ——- by several terms of the expansion in povers of 
x Mat 3 Ma.s ‘ls 


a(x = e[kT), the numerator reduced to a sum of simple Gauss integrals y KE 
e^* dx, whereas the denominator, or, preferably, the expresion in square brackets 


reduces to well- known integrals of the type {re CNET 


On the Mobility uy, 


I xil 
HET 
From Boltzmann's equation 
exp | B n) 
16 2 e? m* kT 
P A m X cu E p NES XIV 
PIE Yu E | «(9 oy : ; €^de | (XIV) 
+ exp Ln 
by eq. (19) 
200 ma... el 
ncs 2'h Agr de (XV) 


l + exp Ge n) 


substituting in (XV) the quantities (XIII) and (XIV) together with t from eq. (X) 


and on appropriate transformations, we have 


pa exp(x— 97 — 
4 J3 
on = $ 2 og a hte Cy E mas’ |l + exp E 2)? 


e (kTy 2 mt vlde (25) 
l + exp (x — 7) 
or 
— hte C 4 
Hoh = 2 Vin ET BT m 5 (26) 
wherein 
| mas!" xh dx 
mu! = 2 l = SED (x 3 n) (27) 


AN E LCR T 
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mas? |l + exp (x — 9)? 
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RECTIFICATION OF PHOTOELECTRIC SPECTROMETER 


By Kazimierz W. OSTROWSKI 
II. Chair of Physics, Mining and Foundry Academy, Cracow 
(Received September 4, 1959) 


In the present paper formulas are derived for the tolerance of the disposition of the 
spectrometer camera exit-slit. Other formulas accounting for the curvature of the spectral 
line set a limit to the admissible height of the exit-slit. A practical rule for the recti- 
fication of the spectrometer is also given, the principal feature of which consists in the 
use of the Hartmann diaphragm. From higher order diffraction, the effect of neighbouring 
lines on the data recorded is assessed. Computation is adapted to the MICII 22 spectro- 


graph. 


1. Introduction 


In direct spectrometry, instead of the photographical plate a photoelectric detec- 
tor (photomultiplier, photosensitive G.—M. counter) is introduced; the latter is 
displaced parallel to the focussing plane, within which displacement of the exit-slit 
occurs. During the measurement, the width of the slit is fixed so as to entirely comprise 
the spectral line. It is but rarely that a slit of lesser width than the line is used. In 
this case continuous recording is obtained with an instrument provided with an 
automatically recording milliampermeter wherein motion of the paper tape is coupled 
to that of the exit-slit with selsyns. 

. Although direct spectrometry has now been in industrial use for about 15 years, 
there exists as yet no systematic treatment of the problem of rectifying the camera 
with the photoelectric detector accounting for such factors as the power of separation, 
curvature of line, ambient temperature, pressure and atmospheric moisture content. 

Some papers (Hanau and Wolfe 1948, Bryon and Nahstoll 1948, Naish 1951, 
1952, Mostyn 1952, Soda 1953, Hagenmah 1954, and Carlsson 1954) give these pro- 
blems passing consideration, and this but from an experimental standpoint. The 
present author intends to complete the foregoing papers by giving computations 
of the admissible tolerance in disposing the camera with the photoelectric detector. 
Such computations may be of interest, as similar problems have been dealt with 


` án the case of grating spectrometers (Minkowski 1942, Rupert 1952, Fastie 1952 a, b.) 
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These authors considered i. a. the problem of chosing the curvatature radii of the 
entrance and exit-slits so as to eliminate aberration and astigmatism at all wavelengths 


when using a long slit. 


2. Assessement of tolerance in disposing exit-slit 


Assume a prismatic spectrometer to have been adjusted as spectrograph. By 
displacing a very narrow exit- slit parallel to the focussing plane, the profile of the line 
can be determined. The latter will be wider then natural for three reasons: 

l. the plane of motion of the slit is at a distance g from the focussing plane, 
which moreover can differ for the various spectral regions, 

2. the exit-slit generally subtends a non- zero angle, y with the wm line, 

3. the spectral line possesses curvature of radius o, 


a. Assessment of admissible distance g between focussing plane and plane of motion 


of slit 


Fig.l shows the situation dealt with. R denotes the line-width measured at g-<0 
assuming the width of the exit-slit, the intrinsic width of the line, and the width 


—————— 
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wherein sin f is a quantity identical with the relative aperture of the camora i. e. 
the ratio of the radius of the objective lense and its focal length. In the present con- 
siderations, values of g will be admitted for which R — d, with d denoting the width 


Fig. 2. Intrinsic width of line, d, and admissible width R as resulting from distance g of slit from focus- 
sing plane 


of the spectral line (Fig. 2). If this width be assumed as the “instrument“ line-width 
i. e. a half-width of the zero order diffraction maximum, then d is given by the well- 
-known formula of the theory of optical instruments: 

gti À sin & . (2) 


sin f 


For the HCII 22 spectrograph used by the present author in investigating photosen- 
sitive G.—M. counters (Ostrowski 1959), the following values of d and &-hold (Pro- 
kofiev 1951): l 


Table 1 
2000 Å d= 5.6 microns & = 47° 
2570 À 9:054 41? 
5600 À 23.08. 5 3% 


According to catalogue, the relative aperature of the objective lense in the MCI 
22 spectrograph is 0.05 for the spectral region of about 2570 À. If these values are 
substituted in eq. (1), an admissible distance of about 0.1 mm is obtained. Let us 
compare this value and that assumed by Prokofiev for the same spectrograph when 
using photographic plates. In this case d is limited by the grain-size (ca. 20 microns). 
Prokofiev assumens the admissible distance of the plate and focussing plane to be 
0.5 mm. It should be stressed that additional widening arises in the plane from dis- 
persion of radiation in the emulsion, especially at greater wavelengths, for which 
gelatine exhibits lesser absorption. Resulting from non-perpendicular incidence, 
the long-wave side of the line will show greater diffluence, Hence adjustment of 
the photoelectric spectrometer requires a greater degree of accuracy; however, its 
separating power is higher, too. Similarly, temperature will affect the precision of 
adjustment to a greater degree (resulting in variations of the position of the line 
and in variations of g). On the other hand, such variations can be set right within 
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a very short time by investigating the profile of the apprioprate line, whereas the 
photographic method makes it necessary to make test photographs. The value of g 
can be modified by introducing appropriataly cut props between the spectrometert 
and the camera; of a set of these are to hand, the camera can be set to the degree 


of accuracy required. 


b. Assessment of maximum admissible angle between slit and line 


Hitherto it has been assumed here that the height of the entrance-slit is so small 
that inclination in the exit-slit has no widening effect upon the of profile the line. Fig. 3 
shows a spectral line L and profiling slit 5. As admissible angle y, the one will be 
here assumed at which the of edge of the slit is the diagonal of the spectral line. Then 


d 
CRAS eia pi (3) 


For d = 10 microns and | = 10 mm, the admissible angle is 0.059. Simulta- 
neous setting of the slit to this degree of accuracy and continual control of the spectral 
line profile is a very tedious affair. The present author proposes to make this task 
easier as follows: an appropriate (narrow) line is subjected to separation with a Hart- 
mann diaphragm placed just before the entrance-slit (Fig. 4). When profiling the 


(e 
D ET 
A aiaphragm 
entrance slit 
Fig. 3. Relative position of line and Fig. 4. Division of line for obtaining sufficiently small 
narrow exit-slit angle y (method of blending of two maxima) 


line two maxima are obtained; these can be made to coincide into one, presenting 
double intensity, if the camera is constructed so as to admit of rotating the slit through 
several degrees of the arc. The method proposed makes observation independent 
of the widening of the line profile arising from curvature. The higher the entrance- 
-slit, the easier it becomes to attain parallel setting of the line and the exit-slit with 
an increasing degree of accuracy. 


3. Limitation in height of entrance-slit relating to line curvature 


The spectral line presents the shape of a parabola whose convexity is directed 
towards the short-wave end of the spectrum. This results form accounting for refrac- 
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tion outside the principal section of the prism. The radius of curvature of the line 
is given by the well- known Kaiser (1900) formula: 
n? f : 
p 2(ni — 1) ctg a (4) 


with f denoting the focal lenght of the collimator, 
n — the refractive index for the line considered, and 
i the angle of incidence, at minimum deviation (in the principal plane). 
Fig. 5 shows the spectral line L, account being taken of the curvature. It will 
be here assumed that the admissible effective height l of the exit-slit is the one at 


a b 
Fig. 5. Relating to definition of admissible effective height of exist-slit (admissible height of spectral line) 


which the tangent A to the first edge of the line becomes the chord of the second 
one. Using a slit whose breadth is double that of the line, it is possible the comprise 
the latter entirely. By the theorem of Pitagoras, 


l = 4yo?— (o— d)? (5) 
Neglecting d? with respect to 2ọd, we have 


1 = AV2- Vod. (6) 

The admissible entrance-slit height will be smaller in proportion to the magni- 

fying factor of the spectrograph (for HCII, this amount to about 1.5). If the exit-slit 

is several times wider than the line, then the admissible effective height is obtained 

by substituting the former's width for d in eq. (6). The paper by Gates (1952) brings 

a detailed computation of g as function of the refractive index, for different angles 
of the prism. The following table has been prepared from Gates' graphs. 


Table 2 >= 
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It is seen that in direct spectrometry the prisms of smaller angle are of greater 
advantage. Computation yields the following results for the height of the exit-slit 
in the case of the MCIL 22 spectrograph: 


Table 3 
o = 35cm d (microns) l (mm) 
f 1 

5 7.4 

10 10.4 

20 14.8 

40 20.8 

80 29.5 

100 33.0 


The present author has carried out the computations for coincidence of the 
optical axis of the spectrometer and the centre of the entrance-slit. Were the axis 
to coincide with the edge of the slit, the admissible height would amount to one half 
of that computed (Fig. 5b). The correction y as proposed by the author automati- 
cally for the inclination of the exit-slit accounts for deviation of the centre of the 
slit from coincidence with the optical axis of the spectrometer. Eq. (6) serves to deter- 
mine experimentally the radius of curvature. For this, we proceed to find the diffe- 
rence in height between two positions of the camera that yields a difference in posi- 
tion of the centres of the profiles equal to the width of profile of the spectral line. 
That difference in height, when read on the scale of the arrangement for vertical 
displacement of the plates, will equal //2. Then eq. (6) will yield the radius of curva- 
ture. From the theoretical point of view the line-width adopted and the way the width 
has beenn defined are not essential. Practically this will affect nothing but the degree 
of accuracy in measuring the radius of curvature. 

Since in direct -reading spectrometry it is of importance that the recorded freq- 
uency of counting be maximum, and since this is in proportion with the area of the 
entranceslit, it seems of interest to consider the problem of an exit-slit with curvature. 
Let o, and og denote the radii of curvature of the slit and the line, respectively. The 
height of xhe slit as in the situation of Fig. 6 will be assumed to be still admissible, 
if the slit and the line have the same width. From eq. (6), with the notation x instead 
of d (the latter is to denote the line-width ontly), we have: 


l = 4y2oyx (7) 
and | 
1 = 4V20 (x + d). (8) 


On eliminating x between both equations, we have 


les d 
l=4 162 : í 
i 01 — 02 2 
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The latter is a generalisation of eq. (6), to which it reduces by substitution of 9, — oo. 
If 9,7205, and denoting the mean value of the radius by o and the difference |0,— oal 


by Ao, we have finally: 
2d 


With d = 10 microns, o = 35 cm and Jo = 2.5 cm (corresponding to the range 
of 2000—2500 À), we obtain / — 3.8 cm. What has been said above concerning eq. 
(6) still holds for eq. (10). 

Finally, let us consider the divergence between a circle and a parabola of the 
same curvature. A spectral line is stricly, a section of a parabole .Tha equations of 


p 
U 
A 
z 
Fig. 6. Relating to definition of admissible Fig. 7. Relating to computation of divergence 
effective height of exit-slit (with curved slit) between abscissae of circle and parabola 


the circle and parabola are of the following form, respectively (cf. Fig. 7): 
yx + Gi — 9)? = 0% 
» —120 Xp. 


The condition y, = y, yields a? + 20(x, — x,) = 0. x, is approximated from eq. (6). 
Finally, the following approximate formula is obtained: 


[4 


MR SITES 


(11) 
With o = 35 cm, the value of the difference begiris to be of significance at l = 25 cm. 
It is seen cleary that the significance of substituting the parabola by a circle of the 
same radium of curvature is negligible. ; 

The construction’ of curved slits is easily carried out by evaporation onto a quartz 
plate and subsequent incision (Polster 1951). The paper by Firestone and Randall 
. (1938) deals with the construction of a curved slit of adjustable curvature. Fastie 
(1953) and Crosswile and Fastie (1954 and 1956) used 10 cm slits of 20 cm radius 
of curvature for gratirig spectrometers. 
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4. Diffractional effect of meighbouring lines on recorded intensity 


It is a well-known fact that a spectral line results from a pattern of diffraction 
on the prism. Diffractions on the slit does not affect the resulting picture of the line 
in a way essential to our. present considerations. The half-width of the zero 
order maximum is termed the “instrument line-width“ and is given by the formula 


cited previously: 


re (2) 
sin f 

wherein @ is the angle of inclination of the beam with respect to the focal plane, and 
sin f — the relative aperture of the camera. The diffractional minima occur at a dis- 

tance of the instrument width. From the Fresnel-Huygens principle, the consecutive 
approximate values of the diffractional maxima can be computed using the formula 


: 2 
- (SE). an 


wherein p = 2zn. The following table contains the results up to order 500. 


Table 4 


TIX 


1.5 9.0x 107* 
2.5 4.3X 1074 
3.5 2.4X 1074 
4.5 1.0x 10-4 
5.5 4.0X 1075 
6.5 1.0x 1075 
1.9 2.5x 10 
8.5 4.0X 1077 


Consider the following example: let the recorded line be P 2136.199 À and the 
perturbing one — the interíse Cu 2135.976 À line. These are wavelengths for which 
the dispersion of the HCII 22 spectrograph amounts to 3.5 A/mm, yielding a sepa- 
ration. of 60 microns for the two lines. The instrument width of 6 microns enters 
60 microns ten times. If the width of the exit-slit is 20 microns, then the slit stret- 
ches over orders from 7 to 12 of Cu 2135.976 À. To compute the number of orders 
comprised, the slit should be projected upon the focussing plane along the beams 
perpendicular to the slit. Hence in the example under consideration not 3 but 5 
orders are comprised. Thus, the effect of the Cu line will amount to something like 
an additional 59/5, of the Cu line intensity. If the last line is of an intensity of e. g. 100 
times that of the P line, then one half of of the effect recorded i is due to diffractiona 
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perturbation. The problem arises of how to account for the effect under consideration. 
. It would seem right to introduce appropriate corrections for the perturbating line, ac- 
counting for its distance from the line measured ,.as well as for other factors. However, 
it would seem more correct to proceed along the following lines; Table 4 shows the 
effect of diffraction to be solely a quadratic diminishing function of the order of the 
spectrum. This effect is, thus, a far-reaching one. A number of neighbouring lines 
affect the line recorded, yielding a certain background. Hence the most correct ap- 
proach consists in giving the effect of diffraction the treatment of a background. 
The background may be considered to be due largely to diffraction. 


5. Conclusions 


` The foregoing cosiderations lead to the conclusions that, although the demand 
for high luminosity in photoelectric spectrometers is consistent, their construction 
should fulfill somewhat different conditions than those observed in the construction 
of spectrographs. Since the slit is to comprise the line with excess, it is possible partly 
to renounce conformily in favour of high liminosity. It may be of advantage to use 
high slits, high 45? prisms (instead of the 60^ prisms commonly in use) having great 
focal lengths, and cylindrical lenses. With these conditions, higher luminosity and 
a greater radius of curvature will be obtained. 
The author wishes to thank Professor Dr L. Jurkiewicz for his critical remarks 
and his valuable discussions. 
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MEASUREMENT OF THE SLOW NEUTRON SPECTRUM OF 
A NEUTRON BEAM FROM THE WWRS REACTOR BY MEANS 
OF A CRYSTAL NEUTRON SPECTROMETER 


Bv D. A. O’Connor J. SOSNOWSKI 
Institute of Nuclear Research, Polish Academy of Sciences, Warsaw 
( Received September 19, 1959) 


The efficiency of a crystal neutron spectrometer was determined experimentally 
using a double crystal arrangement. On the basis of these results the slow neutron spec- 
trum of a neutron béam was determined in the range of neutron wavelength 0.5 — 2.5 À. 
The variation of crystal reflectivity with wavelength was compared with a theoretical 
expression for this dependence. 


The inherently high energy resolution for thermal neutrons of the crystal spec- 
trometer makes it a very convenient instrument for many types of neutron measu- 
rements, such as cross section measurements etc. Difficulties in determining the 
energy dependence of the efficiency of the crystal spectrometer hinder its application 
in precise neutron spectrum measurements and in this respect it is much less con- 
venient than the mechanical chopper devices. Sturm [1] has used theoretically derived 
expressions for the efficiency in a study of a reactor spectrum. In the present work 
an attempt was made to measure the efficiency experimentally. 


Theory 


We consider first a single crystal spectrometer consisting of entrance collimator, 
monochromating crystal, exit collimator and neutron counter. If with the exit colli- 
mator and counter set at an the angle © with respect to primary beam, the crystal 
is rotated to give the maximum intensity of Bragg reflected neutrons, then the coun- 
ting rate of the counter is given by (1): 


P(0) = An(A) 6, cos O f f L (A, à). R(A, f) Ta (A, « — 28) da df (1) 


where n(A)dA is the flux of neutrons entering the first collimator of wavelength A— 
—A+dd; L(A, a) is the efficiency of the first collimator i. e. the probability that 
a neutron of wavelength À entering the collimator at an angle « in the plane of diffrac- 
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tion with respect to the collimator axis, will pass through the collimator; 1,(A, a) 
is an analogous quantity for the second collimator; R(A, p) is the reflection coeffi- 
cient of the crystal i. e. the probability that a neutron of wavelength À will be Bragg 
reflected by a crystal mosaic blok lying at an angle f with respect to the effective centre 
of the mosaic block angular distribution; e; is the efficiency of the neutron counter 
and À is a constant depending on the constant geometrical factors of the system. 
The limits of integration in equation (1) are determined by the geometry of the system. 

In equation (1) only the first order of reflection is included, in general one should 
include higher orders and P(0) will be the sum of a number of similar expressions 
for À, À/2, A/3 etc. 

In the present case it was demonstrated experimentally (see later section on 
results) that 7,(4, x) and L,(4, œ) could be regarded as independent of À and both 
could be represented by the same expression: 

I(a) es 28 

Substituting in equation (1) for /,(&) and J,(2B—«) and integrating with respect 
to æ in the limits —c to oo we obtain: 


P(0) = B.n (4) e, cos O ih R (4, B) e~ Sdp (2) 


where B is new constant which includes a function of Ø. To determine n(A) it is now 
necessary to determine e, and R(A, B) and to eliminate terms due higher order reflec- 
tions. e, can be calculated from the known contants of the counter. 

If the form of the dependence of R(A, B) on B is known, or is known to be inde- 
pendent of 4, then the dependence of R(A, B) on À can be determined experimentally 
using a double crystal spectrometer in the (1, — 1) setting. Consider an arrangement 
consisting of entrance collimator, first crystal, second crystal, exit collimator and counter. 
Neglecting for the moment higher order reflections we.can write for the counting rate 
in the (1, —1) position: 


P, (0) = C. n (å) e, cos O n I, (a) I, (a) da f R2(A, B) df 
= D. n (À) e, cos O f R2(A, B) dB (3) 


where it is assumed that the two crystals are identical. From equations (2) and (3) 
it is seen that if the form of the dependence of R(1, B) on f is known then the varia- 
tion of R(A, B) with A can be found from measurements of the ratio P,, (9)/P(0) 
as a function of 4. In particular if R(A, B) is given by 


R(,B)—R'()e^ 2 (4) 


i. e. the reflection curve has a Gaussian shape, then: 


Pu) 
P(6) 


= constant - R’(A) (5) 
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In fact for an ideal mosaic crystal R(A, f) is given by [2]: 


Pop 
Rag eee (6) 


l+ We 2% 


where y is a measure of the mosaic spread of the crystal and W is a known function 
of À. Using this expression for R(A, B) it is clear that the ratio P,, (@)/P(@) will not 
be a simple function of R(A, f). However the form of equations (2), (3) and (6) sug- 
gests that if for R(A, B) we substitute a simple Gaussian expression of the form given 


by equation (4), the resultant error should be relatively small. It should be noted 
2 
that in equation (2) R(A, f) is multiplied by the Gaussian expression exp — p 


and in equation (3) R(A, B) is raised to the second power which again should give 
a closer approximation to a Gaussian expression. The accuracy of this approximation 
will obviously depend on the parameters of the collimators and crystals used, but 
for the moment we shall continue to use expressions (4) and (5) and return to an 
estimation of the errors involved on the basis of the experimental results in a later 
section. 

When higher orders of reflection are not negligible then the method outlined 
above must be modified. In practice over a considerable range of values of À, from 
0.2 À to 1.5 À, the second order reflection may be neglected. Between 1.5 À and 2.5 À 
the second order component considerably increases but the third and higher orders 
are negligible. Including the second order, equation (2) becomes 


. P(0) = Bn (2) sos 6 | RAP oF did Bn (5) ea cos | a(t P aß 
3 
| = P,(6) + P; (©) (7 


where R,(A, B) and R,(A, f) are the reflecting powers for first and second order reflec- 
tions for a first order wavelength setting of À. The ratios P,(@)/P(Q) and P,(0)/P(0) 
can be determined by measuring the transmission of calibrated 1/v filters (see refe- 
rence 1) for the reflected beam. 

Including the effect of second order contamination, the double crystal measure- 


ment gives: 


Pu(9) = C-n(A) a cos 0 [ n (A, B) dB + Cn (4) Ey/2 COS o fre (A, B) dB (8) 


In fact it was found that the term in RRA, B) was negligible compared with 
that in R? (A, B). Using the value of P,(@) obtained by the filter measurements we 
obtain from (7) and (8): |j 


PHONE DE (9) 
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Apparatus and methods 


Measurements were made on a double crystal neutron spectrometer which has 
already been described [3]. The collimators used were steell Soller slit type of nominal 
half — angle 10 min. Two copper single crystals were used, both having a polished 
surface cut parallel to the (111) plane according to methods which have been described 
[4]. The cut surfaces were of dimensions approximately 100 x 40 mm. Two neutron 
counters were used of identical construction, one filled with boron trifluoride of natural 
boron isotopic composition and the other filled with B1° — enriched boron trifluoride, 
the latter counter being used for measurements at À < 0.5 À. The calculated efficiency 
of these counters is given by 


t l-2 E 1 
£4 =e 3.53 X 10 ZU ETC 6.35 x 10- ay natural boron 


e71.88x10 (] _ e 3.384) 


£4 = enriched boron. 


The source of neutrons was the No 4 horizontal hole of the WWRS reactor. 
The situation of the end of the hole with respect to the core is shown in Figure 1. 
Between the end of the hole and the nearest fuel elements an air — filled aluminium 


reactor core 


fuel element sections 


air cushion 


Moderator (water) beam hole No4 


Fig. 1 


can be inserted into the light water moderator, thus increasing the fast neutron flux 
available from the hole. This air cushion is of effective thickness about 75 mm, and 
when in position leaves a water layer of about 30 mm between the hole and the nearest 
fuel elements. Measurements were made with and without this air cushion, at a reactor 
power of 2 MW and with a moderator water temperature of 25 to 30°C. The stability 
of the reactor power was such that a monitor counter proved unnecessary, in fact 
the day to day reproducibility of results was better than 0.5 7t 

The 1/v filters used were plates of a glass containing 11 96 by weight of B40,. 
Plates of four different thicknesses were available and they were calibrated using 
a singly reflected beam of 1 À neutrons. The dependence of absorption cross section 


was not tested but assumed, which is justfied for À > 1 À and for the boron content 
used. 
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Results 


In figure 2 there are given curves of P(O) as a function of wavelength for the 
range 0.1 to 2.5 À for the two cases (a) with air cushion and (b) without air cushion. 
These curves have been corrected for a background counting rate which was taken 
to be that given when the crystal was rotated through 30 min from the Bragg position. 
The background was relatively constant except for very small Bragg angles and 
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amounted to about 1 % of the peak intensity at the maximum of the curve. The small 
narrow (about 10 min wide) valleys in the curves, produced by parasitic reflections 
in the crystal have been smoothed out. The statistical errors are everywhere very 
small but the reproducibility of the intensities is of the order of 0.5 %. 


Single crystal rocking curves were taken for both crystals at wavelengths of 0.5, 
1.0, 1.5, 2.0 and 2.5 À. The curves obtained are all practically identical as will be 
seen from Figure 3 where points obtained at various wavelengths have been collected 
together on one curve by normalising at the maximum. The half — width of this 
curve is 10.8--0.5 min. This was compared with the curve obteined by rotating the 


- second collimator in the primary beam at very low reactor power, which gave a half- 
width of 18.0+0.5 min. 


Both these curves are well approximated by Gaussian functions and using the 
notation introduced earlier we obtain Ø = 5.64-0.15 min and 7 = 2.3+0.5 min. 
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Fig. 3 
The fact that the rocking curves are practically identical for all values of A justifies our 


assumption that /,(x) and (x) are independent of A. It would further appear that 
the shape of the R(A, B) function is also independent of A, since there is no possibility 
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that changes in the width of the R(A, B) distribution with change in A could be com- 


pensated by opposite changes in the width of the (œ) function, both in fact might be 
expected to increase with À. 


Measurements of the content of higher order reflections in the singly reflected 
beam were made with the boron filters mentioned above in the range 1.0 to 2.3 À. 
Two filter thicknesses were used so that theoretically the percentange of first second 
and third order neutrons could be determined. In fact below 2.0 À the third order 
was negligible. Above 2.0 À the measurements were not accurate enough to give 
the third order with any certainty but the calculated first order intensity was the 
same within +2 % when the third order was assumed to be negligible and when 
it was included. The composition of the beam (with air cushion) as a function of À 
is shown in Figure 4. 

The reflectivity of the crystals as a function of À was measured for several values 
of A in the range 0.25 to 2.5 À. For À < 1.0 À R (A) equation (5) was used as a basis of 
calculation. Above 1 À the results obtained from the filter measurements were 
combined with those from the double crystal measurements using equation (9). 
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After correcting for higher order reflections, and for the total spectrometer 
efficiency we obtain curves for n(A). These are shown in Figure 5 for the two cases 
with and without air cushion together with a Maxwell distribution for an equilibrium 
temperature of 350°K. The portion of the curve below 0.5 A is omitted since for 
shorter wavelengths the accuracy of determination of the efficiency is progressively 
poorer. Both curves show a maximum at 1.044-0.01 A, which corresponds to a Maxwell 
spectrum of temperature 350+7°K. This is about 50°K above the temperature 
of the moderator which is normally in the range 25 — 30°C. Introduction. of the air 
cushion results in an increase in the intensity of faster neutrons and a decrease 
in intensity of the “colder“ neutrons. Above 2.5 A both curves fall below the Maxwell 


spectrum. 
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Discussion of results 


On the basis of the results it is now possible to make some estimate of the pro- 
bable errors inherent in the method of measurement of the crystal reflectivity. 

The fact that the shape of the rocking curve does not change with wavelength 
suggests that the widths of the collimator function and the crystal reflectivity curve 

are constant. 

The absence of the increase with wavelength in collimator function width which 
should result from total reflection of neutrons by the collimator material, is probably 
due to the presence of a thin film of oil on the steel shims of the collimator. 

_ The curve obtained for the crystal reflectivity (figure 6) shows that above about 
0.5 À secondary extinction in the crystal increases rapidly with increase in wavelength 
and the crystal reflectivity reaches a “saturation” level above I.0À wykreślić. This means 
that the shape of the R(A, B) curve must change significantly between 0.5 and 1.0 A. 

However calculations of the R(A, f) curve according to equation (6) and using the 
value 7 = 2.3 minutes, gave curves which, while departing from Gaussian shape 
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considerably for À > 0.5 À had effective widths which corresponded to the exper- 
imental value within the limits of the rather lange, +0.5 minute, experimental error. 
The most important possible source of error in estimating the crystal reflectivity 
is the assumption that P,,(@)/P,(@) can be taken as a measure of the integrated 
reflectitivity of the crystal. The quantity actually measured is given by: 


NIE a seat EAR es (10) 
Put R (å, p) e s dB 


The quantity to be used as the effective crystal efficiency is from equation (2): 


[RG Be er (11) 


Using the theoretical expression for R, (å, B) (6) with n = 2.3 minutes, 
the values of expression (10) and (11) were calculated by numerical integration and 
compared. The dependence of the integrated efficiency (11) on À is given in 
figure 7. It was found that for this particular case the ratio’ of expressions (10) and 
(11) was constant to within 4 96 in the range 0.5 < À < 2.0 A. The simplifications 
inherent in the method would therefore appear to be justified. Comparison of the 
theoretical curve, figure 7, with that actually obtained, figure 6, shows that in fact 
the secondary extinction was considerably stronger than that to be expected theo- 
retically. 

It is possible that the value of used in the calculations is higher than the true 
value. A relatively small decrease in 7 will considerably flatten the theoretical curve. 
The flatter the curve the smaller will be the error resulting from the approximations 
in the method of measurement of the crystal reflectivity, since for large secondary 
extinction the detail of the shape of the R(4, B) curve has a relatively small influence 
^on the integrals involved. Other factors which may possibly explain the marked 
flattening observed are the primary extinction which may not be negligible in a crystal 
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of very small mosaic spread (the mosaic blocks might be expected to be of larger 
dimensions than in crystals of greater mosaic spread) and nuclear absorption, which 
increases with wavelength. 

On the basis of this discussion it is possible to estimate the overall accuracy in 
the measurements of n(A). The possible systematic error due to the method will 
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be a maximum of +2 %. The overall error is calculated to m qure +3 RE a 
0.5 À and 1.0 À and above 1.0 À increases to about + 9 % at 2.0 À the increase 


being due to the possible errors in the estimation of is pose of ovo order 
reflections at the longer Wavelength, 
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COLLISION RECOMBINATION IN TLS 


By J. W. OSTROWSKI 
Institute of Physics, Polish Academy of Sciences, Warsaw 
(Received November 17, 1959) 


The relaxation of photoconductivity in semiconductor microcrystalline TLS layers 
was investigated. The hypothesis of collision recombination of the free current carriers 
in the process of decay of the photoconductivity was put forward and verified. The pheno- 
menological model of recombination: r = «np? was proposed (r — number of acts of recom- 
bination per cm? and sec, & — coefficient of recombination , n, p — concentration of free ele- 
ctrons and holes, respectively). A recombination coefficient x = 10-39—10-9! cm? sec”! 
and a relaxation time of t = 0.5—5 msec were measured. The equation of the recombi- 
nation curve was derived. The temperature dependence of the relaxation time and electric 
conductivity was investigated. The parameters determining the conditions of recombina- 
tion, namely: the concentration of the free current carriers, their mobility, the conduc- 
tivity, the position of the Fermi level and that of the acceptor level, were evaluated. 


1. Introduction 


The hypothesis of collision recombination of the free current carriers in the 
process of decay of the photoconductivity in semiconductor microcrystalline TLS 
layers was proposed and experimentally verified. l 

The collision recombination appears when any free current carrier in a single 
act process passes its entire recombination energy to another free carrier and recom- 
bines either with an opposite free carrier or with the center in the forbidden gap. 

Considering the most simple case of an intrinsic semiconductor wherein the 
concentration of carriers able to recombine is n, the probability for recombination W, 
as for radiative and phonon recombination, is proportional to n, whereas in the col- 
lision mechanism, which is one involving three carriers, we have 


W —n?, 


whence 


Ten, 


As yet, papers on the mechanism of collision recombination are few in number. 
No theory of the mechanism is available, and experimental data are scarce. T. S. 
¥ 


(339) 


Mo. „u prrowekh ^ MEE 
Moss (1953) was the first to recognize the collision model of recombination in order- 
to provide a partial interpretation of his results. 

J. A. Hornbeck and J. R. Haynes (1954) in their paper on the recombination 
of minority carriers with capture centres in p-type silicon, mention i. a. having 
noticed inverse proportionality of the relaxation time and the square of the concen- 
tration of the majority carriers. 

The first to deal exclusively with collision recombination were J. W. Ostrowski 
and L. Sosnowski (1954); the respective measurements were carried out on thallium 
sulphide microcrystalline layers. 

L. Pincherle (1955) mentions a nonpublished notice of 1950 by Fan and Pin- 
cherle where the authors express the opinion that the life-time of the minority car- 
riers may be inversely proportional to the square of the concentration of the majority 
carriers. Pincherle, basing on the Auger effect theory, evaluates the coefficient of 
recombination & at 107% > a > 10-33 cm sec}. 

At the beginning of 1955 appeared a paper (N. Sclar, E. Burstein), dealing with 
collision recombination and the remaining two, possibilities, i. e. radiative and phonon 
mechanisms. The authors carried out their measurements at the temperature of liquid 
helium, at a low concentration of the carriers in the range 109—1.6x 1011 cmr? 

L. Bess (1957) computed the cross section for collision recombination and life- 
-time of the carriers assuming participation of intermediate states situated in the 
energy gap. 

The Proceedings of the Conference on Semiconductors held at Rochester in 1958 
bring the summary by P. T. Landsberg and A. R. Beattie (1958) of an attempt at 
computing theoretically the relaxation time for collision recombination. 

The foregoing review of the literature shows, firstly, that papers, both theoretical 
and experimental, dealing with collision recombination, are as yet few in number, 
and, secondly, that whatever data are available come from widely different. semicon- 
ducting materials such as PbS, Si, Ge, and TLS single crystals and polycrystalline 
layers at various temperatures. 


2. Material Investigated 


The material investigated consisted of microcrystalline, polycrystalline p-type 
TLS layers of thickness ranging from 0.3—1.5 u, evaporated in vacuum at about 
450°C on to the cooled wall of the glass cell in which the sample was subsequently 
enclosed (A. Wolska 1953). 
The TLS layers were investigated with the electron microscope and electron 
diffraction method by Kur and Toruñ. The microscope image revealed the existence 
of two crystalline phases. This was confirmed by the diffraction patterns. Moreover, 
the layers were found to be highly inhomogeneous as to crystal structure. Hence 
it was not possible to give an exact interpretation of the layer structure. The linear 
dimensions of the crystals ranged from 0.06 to 1 = 
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Thus, in the case under consideration, the symbol TLS is more a technical abbre- 
viation than that of a chemical compound. 
The layers exhibit photosensitivity within 0.4—1.2 u, maximum sensitivity 


occurring at 0.83 u. The forbidden gap computed from the long-wave boundary 
of absorption, amounted to 1.1 eV. 


In addition, the layers exhibit a photovoltaic effect of a sensitivity distribution 
resembling that of the photoconductivity. 


The resultant dark resistivity, for various samples, amounted to 0.2—10 MQ. 


nos po — 


S» 8; — 
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3. Phenomenological Considerations 


3.0. Introduction 


Table of Symbols 


concentration of free electrons and holes, respectively, in darkness (at 
thermodynamic equilibrium); 

momentary concentration of free electrons and holes, respectively ; 
steady state concentration of free electrons and holes, respectively, in con- 
ditions of illumination ; 

free electrons concentrations in the presence of background illumination; 
variation in free electron concentration resulting from superposition 
of a light impulse on the background illumination; 

concentration of acceptors; 

concentration of electron-invested acceptors; 

concentration of donors; 

density of states in valence band; 

energy of acceptor level as referred to the top of valence band; 

Fermi level as referred to the top of valence band; 

number of acts of recombination in the dark per cm? and sec; 
number of acts of recombination due to illumination, per cm? and sec; 
total number of acts direct collision. recombination, per cm? and sec; 
probability for collision recombination; 

coefficient of the collision recombination; 

number of free current carriers generáted in dark and by light, respec- 
tively, per cm? and sec; 

relaxation time; | | 

time during which the photoconductivity decays to one half of its maximum 


value (half-relaxation time); 


time during which the additional photoconductivity resulting from super- 
position of a light impulse on the constant illumination decays to one half 
of its initial value; l 

dark and total conductivity, respectively ; 


y 
t 
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9; — conductivity for background illumination; 
x=2 

Co 

0; 
X; = cus 

9o 
d — thickness of layer; 
u — “effective“ mobility of free current carriers; 
J — intensity of light incident on 1 cm? of layer; 
S — photosensitivity, defined as S = do/odl; 
y — number of photons absorbed within range of photosensitivity of layer, 


per energy unit of light; 

light source intensity; 

— distance from: light source to semiconductor layer; 

— the part of radiative.energy of the black-body (of temperature 2200?C) 
that 1s absorbed by the layer wirthin the range of photosensitivity; 

— number of photons emitted per 1 sec by 1 cm? of the surface of the black- 

‘-body at 2200°C; 

energy-radiated per 1 sec from 1 cm? of surface of the black-body at 2200°C: 

— temperature of light source; 


CS 
| 


Su © 
| 


— coefficient accounting for light absorption and scattering along the path 
to the layer; 

time during which the photoconductivity decayed in the measurement 
of the maximum amplitude of photoeffect. 


& 
| 


Let us now approach:the mechanism of collision recombination, confining the 
problem to phenomenological considerations. 

From the facts that: 1)dark conductivity was of the p-type, 2) absorption of 1018 
photons per cm? and sec, the life-time of the carriers being about 10-3sec, raised 
the conductivity to no more than the double of its value in darkness (every photon 
absorbed is assumed to give rise to an electron-hole pair), and 3) the intrinsic carrier 
concentration can amount to about 1019 cm, the following model may be derived: 
In the dark, the 'concentration n of free electrons is very low, whereas that of the 
holes, po, is comparatively high, thus: 


No € po (2) 


Due to illumination, the concentration of the electrons and that of the holes increase 


by an equal amount, which considerably exceeds the value of the electron concentra- 
tion in darkness: 


P— Po =N— Ny =n; 
hence it may be assumed that, in the presence of illumination in the steady state, 


p=n-+ po ~ (3) 


re! 
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Various specialized models for collision recombination that seemed plausible 
with respect to the semiconductor investigated were discussed and tested by the 
experimental data available. Here, account was taken of the possibility of a process 
involving intermediate recombination centres in the forbidden gap. The ensuing 
considerations refer to the model which was found to yield the best agreement with 
the experimental results. 

The case will be considered where the dominant process is one of the recombi- 
nation of a free hole with a free electron, the recombination energy being transferred 
to another free hole in the neighbourhood. The probability for such a collision, per 
unit time and per electron, is proportional to the square of the hole concentration: 


W = ap? E 
The number of acts of recombination per unit time is thus 
r—anp? (5) 


Assume that, in one act of recombination and one act of optical generation, 
one electron-hole pair vanishes and arises, respectively, with every photon absorbed 
generating one pair. 

The equation of the kinetics of the process under consideration is 


d(n + p) 


dt Eo 18; — Sg — oii | (6) 


Here, the notion of “effeetive“ mobility has been introduced on purely practical 
grounds; the meaning thereof will be discussed in detail in Section 3.5; here, let it 
suffice that the relation 


o = u:e(n + p). 


is assumed, justifying the presence of the sum n +p in the derivative of eq. (6). 
The number of photons absorbed per sec by 1 cm? of the layer is equal to that 
of the acts of generation; thus, 


1 
gi —2y J (7) 


in the dark, at thermodynamic equilibrium, we have 


8.” To (8) 


equilibrium of the process is assumed to be preserved in conditions of illuminations. 
By eqs. (3), (5), (7) and (8), eq. (6) of the kinetics assumes the form 


dn 
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The actual concentration of the electrons is equal to the product of the number of 
those arising per sec per cm? and of their life-time: 
1 
n= 2 8; T; (10) 
by eq. (7), this yields 


n=% Jr (11) 


Assuming the “effective“ mobility of the carriers u not to depend on the illumi- 
nation, the following relations are obtained for conductivity in the dark and in the 
presence of illumination, respectively: 

To = eH po (12) 
6 — eu (n+p) (13) 


The relaxation time v, defined as the reciprocal of the probability for the collision 
recombination considered, 


pens 
T 
is given by 
1 
— = ap? 14 
Dp (14) 
The parameter 
» ! 
X= is (15) 


which is directly given by experiment, is a measure of how the semiconductor reacts 
to illumination; assuming constant (at a fixed temperature) mobility of the carriers 
and accounting for eqs. (11), (12) and (13), this may be rewritten as follows: 


Po Po 
The following quantities are directly or almost directly given by experiment, 
and hence are independent of the model adopted: oo, c, d, X, IR Ty» and S. 
With the foregoing quantities, and assuming the model already proposed, the 


quantities Po, 7, 4, &, T will be computed, as well as the following relations which 
can be verified experimentally: 


(16) 


Ty, = fi (Po, X) 
o — fü) a7) 
S =f; (J; X) 


We now proceed to carry out the programme outlined above: - 
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3.1. Relaxation Time 


As we require an expression for the time of decay of the photoconductivity, 


eq. (9) should be integrated at J — 0. On separation of the variables, the latter is 
rewritten as follows: 


mp ^ 
Integration, with the initial condition n(0) = n,, yields 
1 | (n + py n | 1 1 || 
t = — | n 2 — — 18 
apo n (Nm + po) Po Tim + po %TPo 0) 


From eq. (18), an expression for the half-time of decay of the photoconductivity 
T, is readily derived; substituting n(T,,,) — V, n„ and accounting for eq. (16), we 


have 
1 X +3 PX 
Ne nn PR tee ee 
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Eq. (19) admits of experimental verification within proportionality. 

It is noteworthy that, as the light intensity increases, i. e. for greater values 
of X, the time T, decreases. 

From eq. (14), by eq. (16), the relaxation time 7 is obtained as follows: 


4 
X ap(X-r 1)? ; 
Eqs. (19) and (20) yield a relation between the relaxation time t and the half-time 
of decay of the photoconductivity T1, measured experimentally : 


T (20) 


4 
T = Ty, 
h Kane en (21) 
a 2 Eb 
(X — 1) bhri 1X3 


In experimental conditions, the ratio 7 : Ty, differs but insignificantly from unity. 

From eqs. (20) and (21), both the relaxation time t and the half-time v, are seen 
to be inversely proportional to the square of the carrier concentration in the dark, 
in accordance with the theoretical considerations due to Landsberg and Beattie (1959) 

When the semiconductor is subjected to constant, intense illumination — back 
ground illumination — and light impulses are superposed on this background, the 
relaxation time for reverting from the state arising from the additional light impulse 
to the steady state for background illumination will be a function of the intensity 
of the constant background illumination. By analogy to the derivation of eq. (19), 
the relaxation time 7,,4 is now readily obtained. With 


| n(n,) = n; + 15 An | 
and | (22) 
Ninn =n+An 
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and assuming intense background illumination and a weak light impulse: 
n; > poi Po> An, 


we have 


X— Xi X — Xj : 
se gg [o (+) = ae] p 
The half time of decay of the excess conductivity decreases as the background illu- 
mination increases, in analogy to the time T, 
In the presence of very intense background illumination and weak light impulses 
i. e. in conditions when the conductivity in the dark is much smaller than that in the 
presence of illumination (a situation that can be obtained in very sensitive layers 
only) and when the relaxation time t is proportional to the half-time of decay 7,4 
(the decay of the conductivity due to the light impulse is here considered to be ex- 
ponential), a simple, approximate relation t(o,;) can be derived. Namely, by eq. (14), 
with eqs. (3), (12) and (13) (with c; instead of o in the latter one), we have 


t ~ (0; + 00) *, 
whence, with respect to Oo «€ o;, 
T = 067? (24) 


In the paper by Sosnowski and Ostrowski (5), eq. (24) was derived on the ground 
of simplified assumptions; it has been also obtained experimentally (Fig. 7). 


3.2. Conductivity 


An expression for the second relation of eq. (17) will now be derived. 
We now consider the steady state case, i. e. dn/dt = 0. 
Eq. (9) yields J = a(d/y) n(n + po)?, whence, by eqs. (11), (12) and (16), we 


have 


E: «d ps 


ome 


(c? + a? oo — 00, — oj) (25) 


Similarly to the expressions for t, and v, ,, the function o = f, (J) can be verified 
to proportionality. 

For very great light intensities, when it may be assumed that o, € c, eq (25) 
reduces to the simple relation 


Jes (26) 


which is correct for other models of collicion recombination, too. 

This constitutes an easy; experimental criterium for collision recombination. 
There are no grounds for considering this criterium to be reciprocally univocal: 
the present considerations are phonomenological and do not account for all the possi- 
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ble mechanisms of collision recombination. However, the fulfillment of the condition 
of eq. (26) points to the possibility of collision recombination, and, inversely, if col- 
lision recombination occurs, eq. (26) should be fulfilled. 


3.3. Sensitivity 


In the present Section, an expression for the sensitivity 


will be derived. 
Differentiation of eq. (25), on simple transformations and substitution of X — 
= 0/0, yields 
e 8y . 1 
- adp X(X+1) BX—1])’ 


(27) 


which contains measurable quantities nad a constant factor, and, thus, can serve 
as a criterium of the model adopted. 

The sensitivity will now be expressed as a function of the light intensity. Eqs. 
(27) and (25) yield 
qe CHEN 
T X(8X—1) J 
For 1 < X < 2, the first factor in eq. (28) (which depends on X solely) increases 
rather rapidly, whereas for X 722 it varies very slowly, the more so as X increases; 


S (28) 


5 i $7. oh 
hence, for high values of X (in the presence of intense illumination), we have 


(29) 
and, by eq. (26), 
So" (30) 


Eqs. (19), (23), (25) and (30) provide the answer, on an experimental basis, as to 
whether the model adopted accounts adequately for the prosesses occurring in the 


TLS layers 


3.4. Relations between the Principal Quantities 


It is now possible to give some relations between the measurable quantities and 
the quantities po, œ and u. The formulas of Sections 3.0 — 3.3 yield, on elementary 
transformations: 


i J 


(31) 
Petr qup 
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and 
a? (x ee 
“= (oJ DE à > 
"Ce 
2e y Jt (33) 


3.5. Mobility 


It has been hitherto assumed that the mobility of the free carriers does not 
depend on the light intensity, and that it fulfills eq. (13). In polycrystalline layers 
exhibiting also photovoltaic properties, the “effective“ mobility of current carriers 
can. be considered in a double sense: on the one hand, this is the resultant of the 
volume and surface mobilities for intercrystalline transitions, and — on the other — 
the average mobility of the holes and electrons. If volume mobility prevails, we have 


w—=T7*, where VERS 3, 


i. e. a decrease of the mobility as the temperature rises, whereas if barrier mobility 
predominates; the latter increases exponentially with the temperature. It should 
be mentioned here that measurements point to the latter case. 

The problem of how the mobility behaves as illumination varies, remains open; 
on the one hand, the mobility should be expected to vary with the illumination, since 
the intercrystalline barries generally decrease in the presence of illumination; on the 
other hand, the fact that it is the concentrational and not the barrier mechanism of pho- 
toconductivity that predominates in TLS is an indication that the barriers play but 
a secondary role and, hence, that the mobility is independent of the illumination. 
Eqs. (19) and (23) — (30) were derived on the assumption of constant "effective* 
mobility ; hence, if the shapes of the functions 7,,, (X), o(X), S(X) and a, (X;) are found 
to be in agreement with the experimental results, the assumption that the mobility 
does not vary perceptibly in the presence of illumination should also be considered to 
be correct. 


3.6. Temperature Dependence 


It is interesting to know the temperature dependence of the time r and the 
coefficient œ. From the present phenomenological considerations, no conclusions 
can be derived as to the temperature dependence appearing explicitely. The only 
conclusion from eq. (20) is that the temperature dependence of 7,,1s of an implicite 
nature, involving the dependence of the carrier concentration in the dark, py, and 
the parameter X. 

Eqs. (19) and (23) can be rewritten as follows: 


qh. mac seh 
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If the temperature varies, the remaining experimental parameters being kept constant, 
a decrease in the temperature will lead to a rise in X (or X,, too) and, hence, to a dec- 
rease in the values of the functions A(X) and f;(X;); on the other hand, the hole con- 
centration po will decrease, leading to a rise in the factor 1/p5 (the variation in « is 
negligible). Experimental results point to the same order of magnitude for the relative 
variations of the increasing and decreasing factors. Hence it is possible to obtain 
experimental conditions such that if the temperatures are varied within a finite, 
though may be not very wide range, constant values of the time t, and T, , are mea- 
sured. Before a discussion of the results is presented it should be noted that such 
conditions have been obtained in practice. 


4. Measurements 


4.1. Feeding System 


The light source consisted of a 100 W bulb. The light intensity was controlled 
by varying the distance of the source from the semiconducting layer. 


As far as possible, the pulse method was adopted, so that full use could be made 
of electronic devices. 


Rectangular light pulses were obtained by screening off the light with an appro- 
priate mechanical modulator. 


The layer was fed from a DC voltage source consisting of an anode battery; 
the fall in voltage across the layer range from 1 V to 200 V: 


The variable fall in voltage across the oscilloscope input load due to the variation 
in conductivity of the layer illuminated by the flashes from the modulator was, on 
amplification, applied to vertically deflecting oscilloscope plates and brought to 
observation on the screen. In order to preserve the form of the conductivity varia- 
tion, DC current amplifiers were used. The constant component of the fall in voltage 
due to conductivity in the dark was eliminated by appropriate compensation. | 


4.2. The Parameter X 


à 


The parameter X = 0/09 or X; = 0/09 was measured either with a galvano- 
meter connected in series with the layer, by measuring the current through the latter, 
or with an oscilloscope, by measuring the layer resistivity at which compensation 
occured, as observable on the oscilloscope screen by the "zero position of the 
image. 

With respect to the necessity of eliminating long-lifetime processes of the order 
of seconds, the oscilloscope provided the advantage over the galvanometer of reducing 
the time of measurement to tenths of a second. 


D 
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4.3. Half-Times of Decay 


The time t was measured as follows: first, the maximum amplitude of the 
signal on the oscilloscope was measured, the modulator being rotated sufficiently 
slowly for the duration of flasch and darkness to provide for the full development 
and decay of the conductivity variation. Then the duration of darkness was shortened 
so that the conductivity should have time to fall to one half its full amplitude, which 
could be checked by the decrease in amplitude of the image on the screen; the duration 
of the flash, however, had to be kept at a sufficiently high value for the conductivity 
to have time to develop fully, which was obtained by an appropriate shape of the 
modulator. 

The rotation frequency of the modulator was measured by appropriate syn- 
chronisation with the frequency of a pulse generator. Synchronisation was controled 
on the screen of an auxiliary oscilloscope. 

Observation of the decay of photoconductivity with the galvanometer already 
reveals a process consisting in a rather slow decrease of the conductivity lasting up 
to a minute and sometimes amounting to about 10 % of the total photoconductivity. 
This effect was considered to be one of a nature different from that of the recombin- 
ation investigated. The long-lifetime process is eliminated by applying periods of 
darkness of sufficiently short duration 2 at the moment of measuring the maximum 
amplitude of the photoeffect. The values of 9 ranged from 20 to 90 msec. At longer 
times 9, the result obtained ceased to be reproducible; probably perturbation by 
slow processes intervened. 


4.4. Conductivity 


The conductivity in light and darkness was measured by measuring either the 
current through the layer, or its resistivity. The current was measured with a galva- 
nometer connected in series with the layer. The latter's resistivity was measured either 
by the bridge method, or by finding the appropriate resistivity in a compensating 
circuit at "zero^ position of the image on the oscilloscope screen. 


4.5. Sensitivity 


The sensitivity was determined either by graphical differentiation of the con- 
ductivity curve 


o = f (J), 
or by measuring the rise in current resulting from superposition of a constant, ad- 
ditional light pulse on the illumination of intensity J, to which purpose a galvanometer 
connected in series with the layer was used, or else by measuring the height of the 
image on the oscilloscope screen in its dependence on the light intensity. 
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4.6. The Number of Light Generated Carriers 


Computation of the concentration of the carriers of one sign set free by light 
per sec, i. e. of the right hand side of the expression 


l 


amounted to finding all three quantities one by one, with the rather evident assump- 
tion that one photon absorbed generates a single free electron (and a single hole). 

The intensity (in W) of the white light incident on 1 cm? of the layer, accounting 
for absorption and scattering, is computed from the expression 


M 
Hs ink" (29) 
(n = 0.8). 
The number of photons absorbed is computed from the expression 
Z 
y= © (36) 


for T, = 2200°C and the photosensitivity range of 0.4—1.1 u. 

The values of Z, Q and H are to be found in the appropriate Tables. 

The mean thickness d of the layer could be determined to a gross approximation 
only. All direct measurements (as e. g. by optical methods) pointed to considerable 
inhomogenity of the thicknes; this resulted also from pictures obtained with the 
electron microscope. Evaluation of the thickness from measurements of the area 
and mass of the layer yielded results that agreed within about 50 76. The thickness 
of the layers ranged from 0.3 to Lu. 


Fig. l. Schematic arrangement showing a glass cell for temperature dependence measurements. : 


K — tube for flow of liquid (alcohol), V — vacuum chamber for deposition of layer w by 
evaporation, T — thermocouple suspended in tube near layer, E — electrodes. 


ÿ 
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4.7. Obtainment and Measurement of Temperature of the Layer 


In order to obtain various temperatures within the layer and to measure them, 
glass cells of a special shape for enclosing the layer were constructed (Fig. 1). The 
temperature of the layer varied as that of the alcohol circulating in continuous flow 
along the glass wall (support). The thermocouple measured the temperature of the 
alcohol flowing through the tube. 

The temperature of the layer ranged from +-40° to — 60°C. 


4.8. Discussion of Experimental Error 


The highest amount of error was inherent in the measurement of the layer thickness. 
However, because of the inhomogeneous composition of the layer and, especially, 
its polycrystalline structure, no degree of accuracy in determining the thickness 
could essentially reduce the intrinsic impossibility of strictly determining the “geo- 
metry“ for physical processes. 

If the layer thickness is not known with accuracy, none of the quantities referred. 
to unit volume are known exactly. However, their exact knowledge is not of foremost 
importance for the present investigation; what is really important is that their depend- 

ence on the light intensity and the temperature be known as exactly as possible, 
and this, in fact, was measured with a much higher degree of accuracy. 

The principal source of error in all the measurements resided in certain slow 
processes, of which little is known as yet, consisting in a slow rise in the conductivity 
in the presence of illumination and a decrease thereof when the layer is placed in the 
dark. 

The general character of the fundamental relationships remained the same, 
notwithstanding the presence of individual quantitative differences which, as we have 
seen, was not of essential significance. 


5. Results and Discussion 


5.1. Conductivity as a Function of the Illumination 


The measurements of the conductivity in its dependence on the intensity of 
illumination were carried out at various temperatures. These measurements, which 


presented no complications, yielded the first experimental corroboration of the me- 
chanism of collision recombination i 


Since it was the variation of the expression 


f(s) = 0% + 0205 — acl — o? 


that was essential [cf. eq. (25)], this was plotted in arbitrary units on the axis of ordi- 
nates in Fig. 2, whereas the light intensity was plotted on that of the abscissae, also 
in arbitrary units. If the recombination model adopted is correct; straight lines through 


2 Sh 
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Fig. 2. Comparison of f(o) ^ I calculated from eq. 25 (solid lines) with the experimental data (dots). 
Arbitrary units. 


the origin of the coordinates should result. It is seen from Fig. 2 that the experimental 
dots coincide satisfactorily with straight lines. 

In the presence of very intense illumination ip <i, thus yielding conditions 
for which eq. (26) is correct. The results of the respective measurements are shown. 


a 
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in Fig. 3. The curve is clearly seen to tend towards a straight line of slope 1/3, in 


agreement with eq. (26). 


ee IT N. 


L i =! 1 
10020050 200 500 20005000 20000 J 


Fig. 3. Comparison of o ~ I? at strong illumination, assuming i ~ 0, (dashed line) with the experi- 
mental data (solid line). Electrical current i in «A, I light intensity in arbitrary units. 


5.2. The Half-Time of Decay as a Function of the Illumination 


Figs. 4 and 5 bring the measurements of the half-time of decay of the photo- 
conductivity as a function of the parameter X. The continuous line therein represents 
the v, (X) dependence in agreement with eq. (19), representing a factor depending 
explicitely on X; the dots represent direct measurements of 7,,. The theoretical 
curve is determined as to its form only; its position was adjusted to suit that of the 
experimental dots. The agreement is seen to be satisfactory. 

The dependence of the time T, , on the parameter X; was also confirmed by exper- 
imental results. The continuous line in Fig. 6 represents the results obtained in meas- 
uring ty, as a function of X; (the experimental dots are not shown; X; and ty, 
were measured to a rather high degree of accuracy), whereas the dots represents 
the values of the function 


fer ve 
E E ub 1 — ID 


[cf. eq. (23)]; the difference X— X, was determined to a lesser degree of exactitude, 
the respective values being small as compared to those of X and X;. In this case, 
too, agreement may be said to be rather satisfactory, within the limits of experimental 
error. ! 
.. Measurements of T, , as a function of the conductivity o; in the presence of very 
intense, variable background illumination are found to be in rather good agreement 
with eq. (24). The experimental results are shown in Fig. 7; the experimental dots 


‚lie satisfactorily on a straight line of slope —2. It should be noted that eq. (24) is 
an approximation. ; a 
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Figs. 4 and 5. Comparison of ty, (X) from eq. 19 (solid line) with the experimental data (dots), the 
calculated curve being placed on the experimental dots to show the agreement of shape. 
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5.3. Sensitivity 


The sensitivity as a function of the illumination is represented in Fig. 8, the 
ordinates representing the experimental values of S, whereas the function 
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Fig. 6. Comparison of ty, 4 (X ) from eq. 23 with the experimental data (solid line). The den and cros- 
ses represent the values of the function in eq. 23, the difference X—Xi being determined with a great 
error. À à 
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in agreement with eq. (27), has been plotted on the axis of abscissae. If the model 
adopted is correct, straight lines through the origin should result. Here, too, the 
experimental dots are seen to coincide satisfactorily with the respective straight 
lines. Line l in Fig. 8 was obtained at a lower. temperature than line 2 and, hence, 
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50 
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0 10 20 50 700 750 200 250 
fx (x+2)(3x-17)] 7 


Fig. 8. Comparison of S sensitivity (in arbitrary units) vs. [X(X-- 1)3.X—1)]-1 calculated from eq. 27 
(solid lines) with the experimental values for two temperatures (dots). 
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Fig. 9. Comparison of S sensitivity (in arbitrary units) vs. X parameter calculated from eq.27 (solid 
line) with the experimental data (dots). 
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its inclination with respect to the axis of abscissae is greater, in agreement with eq 
(27), as the factor 1/p? now assumes a larger value. Fig. 9 represents the variation 
of the sensitivity as a function of the parameter X; the continuous line corresponds 
to eq. (27) with the assumption that the factor 8y/æp?d is constant; the dots represent 
experimental results. 

For very intense illumination, eq. (28) assumes the form of eq. (29). Fig. 10 
shows the dependence of the sensitivity on the intensity of illumination. 


2 8 20 50 200 500 2000 5000 20000 J 


Fig. 10. Comparison of S ~ 1-1 at strong illumination from eq.29 (dashed line) with the experimental 
data (solid line) in log-log scale. . 


5.4. Shape of the Relaxation Curve 


The foregoing experimental results yield the shape of the relaxation curve for 
photoconductivity : 


o = op f(t). 
From eq. (19), on substituting X — ane zw dj and y — 2n + po and 
Co Po Po 
assuming e. g. X — 4, we have 
yl | 1 3 
t = In 0.6 —— — 2} ————|. 37 
fos y lS (37) 


In order to compare the relaxation curve of eq. (37) with the one for exponential 
decay, an exponential was normalized so as to have two points, namely, the origin and 
the half-time of decay, in common with the curve of eq. (37). The curve is traced 
in Fig. 11. Decay of the photoconductivity is seen to set in at a somewhat higher 
rate than would correspond to the exponential,whereas subsequently it is a much 
slower process. ` | 
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Fig. 11. Comparison of the investigated relaxation curve with the exponential curve. The former first 
falls quicker and then slower than the latter one. 


5.5. Numerical Values of Principal Quantities 


T, = 2200°C 
J —64x10 cm? W 
Z —14x107? 


Q = 2.26X 1021 cm? sect} 
H = 2.12X10? em ? W 

y = xU joule-1 

d =0.5x10% cm. 


With the foregoing numerical values and by eqs. (36), (37) and (38), the number 
of carriers generated by light per sec and cm? of the layer can now be computed at 


Y g; 2x 1018. 
The fundamental eq. (10) yields an electron concentration of 
n,, = 4x 1015 em ? 


for a relaxation time of 2 msec. 
Eq. (31) now yields the darkness concentration po by substitution of the appro- 


priate value of X, at 
po =th6 e106. ema: 


Eq. (32), by substitution of the values obtained, yields numerically the coefficient 


a at room temperature, 


& = 1.4x 10730 cm? sec. 


By direct measurement of the current flowing through the layer and the lat- 
ter's dimensions, the value of the conductivity in the dark at room temperature is 
obtained: 


ay = 5 X 1075 Amp Volt”! cm. 


Eq. (33) yields the value of the current carriers mobility, i. e. the mean "effective" 
mobility of holes and electrons: 


u = 2x 107? cm? Volt sec. 


The sensitivity varied widely according to the illumination and temperature. 
At room temperature, in the dark, with light impulse from a 100 W bulb placed 
at a distance of 100 cm, the sensitivity amounted to 


S = 600 cm? W-1, 


the active surface being 1.2 cm?. 

The results obtained with various samples of different photosensitivity, thickness 
and resistivity in the dark exhibited some amount of dispersion; however, considering 
the diversity of the material used, the latter was, in fact, surprisingly small. 

In Table 1, the approximate values of the various quantities at room temperature 
are assembled. 


Table 1 


1015— 1016 3 


cm” 
Po 1016— 1017 cm? 
[^ 10-30—10-3? cm sec”! 
m 5x 107?—5x 107? cm? Volt-1 sec! 
0, 1071—1075 Amp Volt-! cm! 


ij The value obtained for the coefficient of recombination & lies within the the- 
oretical limits predicted by Pincherle (cf. p. — 346), who evaluated it at 


10722 > a > 10-33, 


Hitherto the literature brings no instance of the coefficient of recombination 
having been determined experimentally. 

The small value of the “effective“ mobility of the carriers is in accordance with 
the hypothesis that the principle obstacle in the path of the carriers moving in the 
electric field resides in the surface effects predominating in polycrystalline material. 
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5.6. Temperature Dependences 


Fig. 12 brings graphs of the temperature dependence of the conductivity, for 
various sample. The latter exhibit but small differences in the inclination of the curves 
-representing the variation of the conductivity in darkness, whereas the curves repre- 
senting the conductivity dependence in the presence of light of various intensities, 
and resembling straight lines, exhibit almost constant inclination. The relations 


x0 SA $82 4230440 IT 


Fig. 12. Conductivity © vs. reciprocal temperature. The curve 1—4 correspond to the dark conductivity 
of the various samples, the curves la — 4a and 1b — 4b to the weaker and stronger background illumi- 
nation respectively. 
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in Fig. 13 are analogous to those of Fig. 12. The linear sections of the o (103/T) depend- 
ences as obtained with four values of the light intensity differing by almost one order 


of magnitude are seen to retain an almost constant steepness. 


mo agaa VEN PT que 


Fig. 13. Conductivity vs. reciprocal temperature. The curve 5 was taken in darkness and the curves F 
_ 5a— 5c at different increasing background illumination intensities 

The tangents of the angles (f) of inclination for the various curves are assembled . 

‘ee ML c a! T " go» S y [a | 7 

in Table 2. The initial section of the c (10%/T) have been considered. x 
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In accordance with the relation oy ~ “po, curves 1 
the temperature dependence of the product 


H a "Po i) (38) 


whereas curves 1 a, b — 5 a, b, c in their final, linear part illustrate that of the product 


l i! 
TEN ^ 


n > Po: 


5 (Figs. 12 and 13) illustrate 


since now we have 


which leads to o ~ un, and with respect to n~ T ((10), (11)) we have finally eq. (39). 
which holds for a given sample and a given intensity of illumination. 

From eqs. (38) and (39) it is seen that the functions po(T) and p(T) are obtained 
if the function t(T) is known experimentally. 

Fig. 14 illustrates the temperature dependence of the relaxation time v. The 
upper curve resulted at a slightly larger value of X i. e. a somewhat more intense 
light impulse, than the lower one, a circumstance of no great importance. The relax- 
ation time is seen to descrease linearly as the temperature rises. 
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Fig. 14. Relaxation time Tt vs. temperature. The upper curve corresponds to the stronger pulse illum- 
ination. y 


Since the variation of t(T) with the temperature is linear and rather slow, the 
character of the temperature dependence u (1/T) of the mobility ensures the ex- 
ponential shape of the product in eq. (39). Fig. 15 illustrates the experimental tem- 
perature dependence of the mobility on the basis of the results yielding eq. (39) 
and the dependence of the time T on the reciprocal value of the temperature. The 
mobility is seen to fall exponentially with the reciprocal temperature, a fact that — 
together with the low numerical value — may be considered to confirm the hypo- 


y 
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thesis of the predominant effect of the surface barriers at the microcrystalline bound- 
aries. 

The form of eq. (38), together with the exponential decrease of the mobility with 
1/T (Fig. 15), yields the temperature dependence of the concentration in darkness pg. 
Restricting to the sample of Fig. 13, the po(1/T) function is found as follows: 


1.91 — 0.56) x 10° 3.1 x 10° 
Bells me care 2) 
wherein 1.91 is the angular coefficient of the strainght line in 5 Fig. 13, and 0.56 — 
that of u (103/T) in Fig. 15. From the shape of pọ (10%/T), the position of the acceptor 


states accounting for the conductivity in the dark can be evaluated. 


igu 
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Fig. 15. Mobility u vs. reciprocal temperature calculated from eq. 39 and from t (T) measurements shown 
in fig. 14. 


Let us now consider the problem from a more general point of view. The holes 
in the valence band are assumed to arise from electron capture by the acceptors; 
the latter; moreover, undergo investing by electrons from donors, although these 
are present in a number far inferior to that of the acceptors. 

The corresponding formulas are: 


n, = N, exp = X. E ) (41) 


Po = Np exp E e (42) 


Let us consider two cases: First assume all acceptors to be invested by electrons 
from the valence band, n, = po; thus, there are very few donors, so that 


Po > Ny. i | (43) 
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With n, = py and by eqs. (41) and (42), we have 


whence 
thus, eq. (42) yields 


VEN, N, E, 
po = V2Np Na exp 5,7 (44) 

2kT 
Now assume the other case, namely 

Po € Na (45) 
here, 
n,*spe Nac Np 
whence, by eq. (41), 


E; F 
2N, exp = au) . exp p — Na, 
and eq. (42) yields 


Na E, 
Po = 2Np ND g 23 (46) 


From eqs. (44) and (46), the acceptor activation energy for the first case is seen 
to amount to the double of that for the second. There is no possibility of deciding 
between the assumptions of eqs. (43) and (45) by direct experiment since neither 
the concentration N, of the donors nor that of the acceptors N, could be made acces- 
sible to measurement. 

We shall now attempt to prove that case I does not occure. For this, the right 
hand side of eq. (44) will be computed: 


. "a 
N, = Aer = 2.38x102; 


in the foregoing, the effective mass of a hole is assumed to be equal to the mass of 
a free electron, and T = 290°K. 3 


PE s 0.54 17° a 
SP (- | EP | = s] ml e 
whence, with po = 5x1016, we have 


N, = 10°? enr 3, 


a result that cannot be retained in view of the fact that the number of atoms per cm? 
of TLS is evaluated at the order of 10??. 
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Let us now consider case II. Eq. (46) yields 


Na « 


a number of acceptors 35 times greater than that of the donors, which is quite rea- 
sonable. With respect to the condition of eq. (45), it seems permissible to put e. g. 
N,= 3x10" and N, = 10”, these being figures that are within reasonable limits. 

Thus, eq. (46) is the one to be retained; it yields the following evaluation for 
the position of the acceptors: = 


E, = 0.21 eV; 


this is an approximation, as resulting from earlier evaluations (as, e. g. of po), and 
additionally because of the lack of information as to the value of the effective mass 
of the holes in the valence band. | 

By eq. (42), assuming the resultant; effective mass of the holes and electron 
to be equal to the mass of a free electron, the position of the Fermi level F et room 
temperature can be assessed as follows: 


Foz 0:156 VR à 
If a hydrogen type model of the acceptors is assumed, the dielectric constant 
of the material investigated can be evaluated; this leads to = x 
13.5 tà dd 
r C= 027 = te " a 


which is a reasonable value. l 
Thus, the following approximate model results for darkness.and room tempe- 
rature: "o 
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5.7. Conclusion 


In Table 3, numerical data obtained with the layer with which the greatest amount 
of measurements had been carried out, and which could be considered to present 
average properties are assembled. The concentration n, and the relaxation time T 
are given for illumination with a 100 W bulb from a distance of 100 cm. 


Table 3 


4x 1015 


1.6x 1016 


approximate values cm? Volt-! sec”! 


Amp. Volt?! cm-t 


Values assessed on the 
assumption of m* = m, 


Reasonable values 


3x 1017 


The model of collision recombination as adopted is, in general, confirmed by 
experiment. Nevertheless, it is certainly no more than an approximate one. It seems 
that material of the kind used is hardly adapted to the verification of the finer, quant- 
itative consequences resulting from more detailed assumptions in the model. 

In semiconductor material of this type of structure, having a vast surface and 
an inhomogeneous composition with many impurities, a number of local energy 
levels situated at various depths of the energy gap should be expected to exist. It is 
thus hardly feasible to eliminate entirely the possibility of recombination of the 
carriers through the centres. However, the experimental results seem to indicate that 
recombination through intermediate states is not predominant here. May be the time 
of intermediate recombination is longer than that for the direct process. The long- 
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lifetime processes observed may well consist in the investing of traps of some 
kind or other. 

Further research work on the mechanism of collision recombination should 
be conducted with material of a type free from the disadvantages of the polycrystal- 
line one investigated. In the first place, monocrystals of high photosensitivity, a high 
concentration of free carriers and a well-determined distribution of impurity levels 
should be recommended. 
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DIELECTRIC PERMITTIVITY VARIATIONS IN BaTiO, SINGLE 
CRYSTALS AND CERAMICS RESULTING FROM HYDROSTATIC 
PRESSURE: 


By J. Kurmowski AND J. PIETRZAK 
Institute of Experimental Physics, A. Mickiewicz University, Poznan 
(Received December 7, 1959) 


The effect of hydrostatic pressure on the dielectric properties of BaTiO; monocrys- 
tals and ceramic below and above the Curie point was investigated. Curves of the per- 
mittivity versus the temperature, for various voules of the pressure, were obtained. The 
maximum hydrostatic pressure was 2.000 kG/cm?. The reciprocal value of the dielectric 
permittivity, 1/e, was found to decrease linearly below the Curie point and to rise linearly 
above the Curie point, as the pressure increased, both in monocrystals and in ceramic. 
The Curie point and Curie-Weiss temperature (i. e. the temperarute of the catastrophe for 
the cubic phase) decrease linearly as the pressure rises, at a rate of —4.1x 107$ °C/at and 
—3.6x 10-3 °C/at for the crystal and —4.5x, 1073 °C/at for ceramic. The Curie-Weiss 
constant is not modified by pressure. The maximum value of the permittivity of 
a monocrystal rises with the pressure and the &(T) curve becomes steeper; in ceramic, 
e decreases and the curve flattens down. The results obtained within the region above 
the Curie point were compared with Devonshire's theory. Using the authors’ experimental 
results, the coefficient of voluem electrostriction 7 was computed. Moreover, dalayed 
effects resulting from hydrostatic pressure pulses were investigated. 


Introduction " 


A number of papers dealing with the effect of hydrostatic pressure on ferroelectric 
properties have been published. Vul and Vereshtchagin (1945) found that the capacity - 
of a polycrystalline BaTiO, plate increased as a result of hydrostatic pressure and 
the relative variation of the capacity per unit pressure amounted to 1.2 x 10-5 cm?/kG. 
Koch (1958) investigated time variations of the dielectric permittivity in a BaTiO; — 
single multi-domain crystal and ceramics by pulse-wise applying and removing 
of hydrostatic pressure. Merz (1950) found that the Curie point decreased linearly 
by —5.8x 107? °C/at in a monocrystal subjected to hydrostatic pressure. 

The present investigation deals with the effect of hydrostatic pressure on the 


1) The present investigation was subventioned in part by the Institute of Physics of the Polish 
Academy of Sciences. 
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dielectric permittivity of BaTiO, monocrystals and ceramics both below and above 
the Curie point. The results obtained above the Curie point are compared with Devon- 
shire’s theory. 

1. Experimental Method 


a) Preparation of Samples 

BaTiO, monocrystals and ceramic were used for investigation. The crystals 
were grown by Remeika’s method, and had the shape of plates 0.2 — 0.3 mm thick 
of an area of 4 — 10 mm?. Silver electrodes were formed on the surfaces by evapo- 
ration in vacuo. Only monocrystals presenting a rectangular hysteresis loop and a die- 
lectric permittivity of about 300 at room temperature were used in the investigation. 

Ceramic samples were prepared from polycrystalline barium titanate by well- 
-known methods (Trzebiatowski and co-workers, 1952). Podwered BaTiO, was com- 
pressed at 4 tons/cm? to form discs 6 mm in diameter. The latter were sintered for 3 
hours at 1380°C, yielding material of low porosity. On polishing the discs down to 
a thickness of 0.6 — 2 mm, silver electrodes were chemically formed on the surfaces. 


b) High Pressure Apparatus 

The apparatus for producing high pressures is shown in Fig. 1. Hydrostatic 
pressures of up to 2.000 at were produced within a piezostat by a screw compressor 
and 1:20 pressure multiplicator filled with paraffine oil. 


Fig. 1. Apparatus for producing hydrostatic pressure, with piezostat 


és u Ze 
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The pressure within the piezostat was measured with a tube manometer, the 
aegnracy amounting to +10 kG/cm?. The electric leads connecting the crystal or cera- 
mic sample with the measuring device passed through the plug in the piezostat (Fig. 2) 
giving rise to a capacity of about 7 pF. The temperature was measured with a thermo- 
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Fig. 2. Plug with electric leads: 1. plug, 2- leads, 3. insulation, 4. sample, 5. insulation (mica) 


couple situated at a distance of about 1 mm from the sample on the side of the groun- 
ded electrode. The thermo- electric force was measured with a Diesselhorst compen- 


sator. 


c) Electric Measurements 

The dielectric permittivity of BaTiO, monocrystals and ceramic was measured 
on a Schering type bridge. The device is shown schematically in Fig. 3. The accuracy 
in measuring the capacity amounted to +0.2 pF. The measuring field did not exceed 
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10 V/cm, and the measuring signal had a frequency of 2 kc/sec. Each measurement 
of the capacity versus the temperature at constant pressure Was carried out by cooling 
the sample from 180°C. The rate of cooling amounted to 30°C per hour and was 
constant throughout each measurement. The variations of the stray capacity resulting 
from changes in temperature and pressure did not exceed +0.1 % and were practically 
negligible. 


AMPLIFIER 


RC 


GENERATOR 


to piezostat 


Fig. 3 Scheme for measuring capacity 


2. Results 


A. Fig: 4 shows the dielectric permittivity of a monocrystal as plotted versus 
the temperature; at various fixed values of the hydrostatic pressure. As the pressure 
rises, the Curie point (the temperature corresponding to maximum e) shifts linearly 
towards lower temperatures by —4.1X 103° C. cm?/kG, and the maximum value 
of the permittivity rises monotonically. The Curie-Weiss temperature © at pressure 
o (the extrapolated temperature at which eis 0) decreases as the pressure rises; 
however, the decrease is smaller than that of the Curie point, and amounts to —3.6 X 
x 1079 °C : cm?/kG (Fig. 5). The Curie-Weiss constant, which amoumts 1.1 x 10* °C, 
undergoes no modification. Hence, as the hydrostatic pressure rises, the interval 
between the Curie point T; and the Curie-Weiss temperature © diminishes. 

Measurements of the dielectric permittivity as a function of the temperature 
and hydrostatic pressure carried out on BaTiO, ceramics (Fig. 6) proved a rise 
in the pressure to produce a linear decrease in the Curie point and Curie-Weiss tem- 


x 
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Fig. 4. Permittivity of BaTiO, crystal versus the temperature, at various values of the hydrostatic pres- 
1 sure: 


1 0 kG/cm?, 
2— 500 kG/cm?, 
3 — 1.000 kG/em?, 
4 — 1.500 kG/cm?, 
5 — 2.000 kG/cm?, 
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5000 


4000 


3000 


1000 


Fig. 6. Permittivity of ceramic versus the temperature, at various values of the hydrostatic pressure: 
ET 0 kG/cm?, 
2— 740 kG/cm?, 


3 — 1.600 kG/cm?, 
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perature amounting to —4.9X 10-3 °C + cm2/kG. (Fig. 7). This value exceeds the one 
for a crystal and is in agreement with results obtained by Jaffe and co-workers (1957). 
The Curie-Weiss constant is pressure-independent and equals 0.9 x 104 °C. The 
essential difference between the results obtained with BaTiO, monocrystals and 
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Fig. 7. Temperature dependance of l/e for ceramic, at the hydrostatic pressure values of Fig 6. 


with ceramic manifests itself in the region of the Curie point: at higher hydrostatic 
pressures, the &(T) curve for ceramic exhibits a marked flattening at the transition 
with a lowering of the dielectric premittivity. 

The authors propose to explain the effect by the fact that ceramics are inhomo- 
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geneous systems consisting of crystallites; internal anisotropic stresses arise in the 
latter when hydrostatic pressure is applied. 

In Fig. 8, the reciprocal value of the permittivity, 1/e, is plotted versus the hydro- 
static pressure at fixed temperatures, for a monocrystal and ceramic. The reciprocal 
l/e is seen to be a linear function of the hydrostatic pressure. Below the Curie point, 
1/e diminishes linearly as the pressure rises, the steepness being temperature depen- 


110° 


ce) 


1120 
E . 


T-150°C 


1100 


60 
i E 
"m 


cM 800 EE 1600 —— à KGcm) © 400 800 1200 1600 —— à (kG/cm?) 


Fig. 8. Dependance of 1/e on the hydrostatic pressure. 


dent, whereas above the Curie point 1/e increases linearly with the pressure and the 
steepness is constant and temperature- -mdependent. 

B. Pulse hydrostatic pressure was applied to and removed from BaTiO, single 
(c — domain) crystals and ceramics below and above the Curie point, the time varia- 
tions of the capacity and temperature of the samples being measured simultaneously. 
Relaxation effects involving the dielectric permittivity were found to occur only 
in ceramic samples below the Curie point. Above the Curie point; variation in time 
of the dielectric permittivity of BaTiO, ceramic samples results only Gom the variation 
| in time of the temperature. 
| When hydrostatic pressure is suddenly applied to the sample or suddenly removed 
therefrom (Fig. 9), the permittivity rises abruptly and, subsequently, decreases 
in time. At 20°C the variations of the permittivity are greater and require longer 
time for decay than at 70°C. ! 

The effect observed greatly resembles the delay effects in an electric field below 
the Curie point as observed by Partington and co-workers (1949) and by Piekara 


. and Pająk (1952). : 
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Fig. 9. Delayed variation of relative permittivity Ae/e in a ceramic sample, as computed without accoun- 
ting for changes in temperature, a) on applying suddenly, b) on removing suddenly hydrostatic pressure. 
Dashed line shows time variation of the temperature as arising from sudden change in pressure. 
graph 1 — at 20?C, 
graph 2 — at 70°C. 


The authors propose to explain the delayed effects resulting from pressure or 
from an electric field by the complex domain structure of the ceramics (presenting 
domains with 90° and 180? walls). 


3. Comparison of Theoretical and Experimental Results 


Devonshire (1949) expressed the free energy of barium titanate as a function 
of the polarisation and strain. For the non-polarized cubic phase, the function is of 
the following form: | 


Aix Bus) = ze (x + yy + zi) Cia (Vy Ze + 22 £s + ox. Yy) + 
+ ea (+ Ys + 2) tix” (+ P2 P?) + 
(PI. Pj PD +3 Ha (PSP? + PP? + pip) + 
+ En (s P yy Py + ze PP) + sins (P5 + P2) +y,(P2+ P3) + 
+ zs (P + P5) + ga (Ye Py Pe + 2x Pa Po + xy Py Py) + 
+ 30" (Pi + PS + PS) Q) 


ATC | 


mr 
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Let us now consider the free energy as dependent on the quantities x,, yy, z, and P, 
alone (for brevity, strain and polarization will be denoted by x, y, z and P, respectively) 
The remaining components of strain x,, y,, 2, and the polarisation components IE, 
and P, will be assumed to vanish; this yields (Forsbergh, 1954) 


A (x... P) = beg (x? + y? + 2%) + cra (yz +22 xy) + By z PPT 


+ Ep (ry) PR" PIE" PAS LPS @) 
It will be assumed that positive stresses correspond to compression, and nagative 
ones — to stretching. Thr stress components are the following: 
OA © A 
X= a d (3) 
9x dy Oz 


For hydrostatic pressure, eq. (2) yields the following expressions: 


0 eX es y) gr P^ 
— 6 = cu y + ey (2+) + gs P? (4) 
— o = n Z + oa (y +2) + e P’, 


c denoting the pressure. The system is solved in x and z assuming that in an electric 
field applied along the z-direction there is transition of the cubic into the tetragonal 
phase; thus, the x and y strains are equal, and the system (4) reduces to two equations. 
The solutions are the following: 


d Cia £12 — Ci? Eii 
x-—y--— o - = Pr (5) 
^ € + 261 (ei — C2) (en + 2032) 
p 1 2€19 812 — (Cu + C13) £11 p? (6) 


— ——;— 6 

C + 2612 (c3 — C12) (Cra + 26 
The free energy A (X, ..., P) can now be expressed as a function of the stresses and 
the polarisation by applying (Tisza, 1951) Legendre's transformation: 


P Ae EA 
VUES cape Cou E 2 See (7) 


Substituting (5) and (6) in eq. (7) and ordering by growing powers of P yields 
A(o,P)=— xot + no Pt ay PE EEG Pra Ph 


or 
A (0, P) = — 4 xo? + no P? + A (0, P) (8) 
. bili _ _ Su +28» 
wherein x — GNOME is the coefficient of compressibility, 5 = ences 


— the coefficient of volume electrostriction, and A(0, P) — the free energy expressed 
as a function of the polarisation at zero stresses. 


y 
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From eq. (8), the dielectric permittivity can be computed: 


934(o,P) 4m _ à. 0*AQ. P) 
Apr C INCL MS DI 
whence 
Jr g + A (9) 
CRT € 


wherein e, is the permittivity under hydrostatic pressure, and € — that at zero stresses. 
By eq. (9), it is seen that above Curie point 1/e is a linear function of the hydrostatic 
pressure. This result is in agreement with experiment (cf. Fig. 8). Tt is now possible 


to compute 9 numerically, assuming = — E — 5,3 x 10-5 (from the diagram of 


Fig. 5) at o — 2.000 kG/cm?, which yields a value of 
dise 175910 35: em? dym 


On the other hand, 7 can also be computed from the relative variation of the volume 
and the variation of the spontaneous polarisation at the Curie point. Differentiating 


eq. (8) with respect to the perssure and assuming o = 0 yields = = 4) RA whence; 
with ES = 49x 10~ (Kay and Vousden, 1949) and P = 18 uC/cm? (Merz, 1953), 
a value of n = 1.7 x 10-1? cm?/dyn is obtained. The relation between the shift in the 
Curie-Weiss temperature and the hydrostatic pressure is obtained from eq. (9) and 
the hydrostatic pressure is obtained from eq. (9) and the Curie-Weiss law: 


0, = 9+ko (10) 
wherein ©, is the Curie-Weiss temperature at pressure o; © — the Curie-Weiss 
temperature at ø = 0; C — the Curie-Weiss constant, and k = —2nC. Substituting 
herein the values of 7 and C from the authors, measurements, a value of k — —3.6x 


x 1073 °C : cm?/kG is obtained, in good agreement with direct measurements of the 
shift in the Curie-Weiss temperature produced by hydrostatic pressure, which yield 
k = —3.6 x 10-3 °C + cm?/kG. 


4. Conclusion 


| Both in monocrystals and in ceramic, the reciprocal of the dielectric permit- 
tivity diminishes linearly below the Curie point and increases linearly above the Curie 
point with increasing hydrostatic pressure. 

The Curie point and Curie Weiss temperature decrease linearly as the hydrostatic 


pressure increases; the steepness amounts to — 4.1 x 10-3 °C/at and — 3.6 x 10-3 °C/at i 


~ 
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respectively, for a monocrystal, and to —4.5 x 10-? °C/at in both cases for ceramic. 
The Curie-Weiss constant does not depend on the hydrostatic pressurre. 

The peak of the dielectric permittivity at transition temperature rises with incre- 
asing hydrostatic pressure in the case of a crystal, whereas it diminishes and becomes 
diffluent in that of ceramic. 

The results obtained above the Curie point are accounted for by Devonshire's 
phenomenological theory. 

Sudden applying or removing of hydrostatic pressure produces delayed variation 
of the dielectric permittivity in BaTiO, ceramics only below the Curie point; a single 
c-domain crystal presents no such effect. 

The authors wish to express their thanks to Professor Dr A. Piekara for directing 
the present investigation and for his constant advice and valuable discussions. 
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SOME REMARKS ON COLLECTIVE DESCRIPTION OF 
ELECTRON INTERACTIONS. THE PARTITION FUNCTION 


Bv JAN SLEDZIK 
Institute of Physics of the Polish Academy of Sciences, Department of Ferromagnetics, Poznan. 
(Received December 8, 1959) 


The nearly-Hermitian variant of collective description is formulated. The correct 
partition function is obtained. Thereby the exact proof is given that the subsidiary con- 
ditions proposed by Kanazawa lead to false results, whereas those of Bohm and Pines 
yield the exact statistical function, provided the renormalizing procedure proposed 
by the author in a former paper is applied to their theory. 


1. Introduction 


In a former paper (1958a) I have formulated à renormalizing method leading to 
elimination of divergences from the collective description. This has been achieved 
by a reformulation of the subsidiary conditions. They differ from the form given 
by Bohm and Pines (1953b) by the additional term iw ,(1—%,)/e(1 +4) (div A(r)), 
where x; are renormalization parameters and their values are restrained to the interval 
0 x x, < 1, with w, denoting the frequency of collective oscilations, € — the speed 
of light, A(r) — the longitudinal vector potential and 4 — the Fourier expansion vector 
(reciprocal lattice vector). 

The difficulties of the collective description are the same as in quantum field 
theory. In either case the state vectors (in the latter case for free fields) are not nor- 
malizable. Hence, in either case, renormalizing procedures are required. Reno- 
rmalizing of the state vectors of the Bohm and Pines theory includes transition 
with dimensionless para- meters %4 to their values near unity. 

In one respect the collective description differs from the picture of two strong 
interacting fields in quantum theory. In facts we know the exact eigenfunction both 
of the Bohm and Pines Hamiltonian and of the subsidiary conditions. Such eigen- 
vectors depend on x;. When x, is nearly unity (strongly coupled plasmon and electron 
fields, because of the divergence of the averaged interaction term (the plasmon part 
represents the infinite superposition of the boson creation operator pairs (a5 a5) 

| , (383) 
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here s is a natural number, whereas with x, = 0 (weak coupling of electron and 
plasmon fields) the plasmon state vector is in reality the vacuum vector. This 
problem will be discussed subsequently and it will be pointed out that the form 
of the Kanazawa theory (1957) should be correct at high temperatures, that is when 
the linear formalism of the collective description ceases to be useful. Hence, it seems 
that the Kanazawa's subsidiary conditions are not correct. 

The physical sense of the ground state function obtained in (1958a) is not the 
same as in Bohm and Pines paper (1953b). This results from the fact that Bohm 
and Pines use a function that does not satisfy the subsidiary conditions. Notwith- 
stanting this, the ground state vector derived in (1958a) is physically correct. The 
reason for it is the following: in quantum field theory the ground state is the state 
of vacuum, provided the interaction of the fields is weak; in the case of strong coup- 
ling we do not know the exact state vector, hence we are unable to say anything what- 


sover about the ground state. On the contrary in collective description, notwith- 


standing strong interaction, this state is obtained with precision. The state thus obtained 
is indeed the ground state, although it differs from the common picture of the ground 
state, i. e. the vacuum vector. 


2. Nearly-Hermitian collective description 


Let us write the extended Hamiltonian and the subsidiary conditions according 
to Bohm and Pines (1953b): 


À aal 4r |^ : 
Hoa = H + 2 ps DkP-k — EJ Ok Dk; (2.1) 
IE| « b, KI<k, 
- 32/2 27 Ok 0-4 — ne? 
=} Bim + ITIN aie? (2.2) 
me FER) 
pO —0, VAE, (2.3) 


wherein q,, p, are coordinates and momenta of the long-range Coulomb field, o, — 
the Fourier transform of the charge density,e — the electron charge, n — the number 
of electrons and V — the box volume. 


The quantities q,, p, and 0, satisfy the relations 


Qi = cx» AFP 0: (2.4) 


the longitudinal vector potential is 


> . f4xc?\ ^ ghee 
án =i | 7 ) Yaz CHAR | (2.5) 


Iki<k, 
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and, respectively 


= I An^ Tee % TE 
jn) 3 3; A0) =—1 E 2 D. eke =i i pas eit T, 
IKI<k. IEI<k. 
(2.6) 
(The transition q, = p-p requires exact investigations. It is admissible when the 
Hamiltonian is stricly oscillatory, or if the electron — plasmon interaction term. 
is small, which is not the case in the Bohm — Pines theory). 


It can be readily shown that the vector potential is not univocally defined. In 
order to prove this let us write the oscillator Hamiltonian in the form 


Hoss Dre Ok dy 1-4): (2.7) 
EI<k, 
The operator (2.7) in invariant with respect to real and orthogonal rotations in the 
plane py, € Q-&^ 


D, —CoS 4} Duet: w, Sin. X q-k> 
5 I ‘ 
Jr cR sin &k Pr + COS V Q-k» (2.8) 
k 
where æ, denotes the angle of rotation. Inversion of the equations (2.8) yields 
Dass 00" Op py — Op SN Ex TR 
rl 29 
dan DUI a Sarah (2.9) 
k 


hence, with respect to the relations 


| i A 
qi 1 Pa (ax — ar) (2.10) 
the subsidiary conditions (2.3) take the form 


— 


Ya 
(22) [cos Ok (a-k + ax) — i sin &k (a-k = ak)] ð= 0, Li < ke (2.11a) 


or 


(er a.., + ek a) B = 0, - n eux (2.11b) 
The equation (2.11b) is easily solved yielding 
o = II D, z Fé (ri; DT Ta) (2.12a) 


[kl s 
k:>0 
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wherein 
Ó, e—(cos 2a,,+i sin 2ay)a ar, D, à (2.12b) 


Hy (ty 0%) = Ex (ty, 5%): (2.120) 


The operator H is given by (2.2); a,, a, are the oscillator destruction and creation 
operators, w, — the oscillator frequency. The vacuum vector is defined by formulas 


$:0,—1, a,0,—0, Dia; =0. (2.13) 


By (2.12b) the orthogonal rotations do not affect the norm of the state vector, though 
they modify the form of the vector potential. Thereby the situation resembles gauge 
invariance in quantum field theory. 

It is known, Kanazawa (1957), that the subsidiary conditions (2.3) or, with respect 
to (2.10), 


(a, a) 0,20, (RE, (2.14) 
are satisfied by the non-normalizable vectors 
D, = ek rd, (2.15) 


This difficulty can be avoided by a renormalizing method entirely different from 
that given in (1958a). For this purpose we retain the original form of the Bohm and 
Pines theory, but define the operators g, and p,in a somewhat different manner. 
We write 


2 ^ 7 An dm v 
Hen =H + 5 2 P ba M ES Ok Pr» (2.16) 
lAI<k, [k| <ko 
Pr D = 0, D (2:17) 
A ho % ^ Ak 
Êr = EJ (a-r + ài). (2.18) 
EN UN. i 
ai (à ) (a), (2.19) 
wherein 
^ I Ak * * 
erg: ak, Ge = a (1 — vk ak a4), 


=’, LL (2.20) 


The newly defined operators are no longer Hermitian. The degree of their deviation 
from the Hermitian form depends on the number VAs v, — 0 the deviation vanishes. 
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We have 
EM ee EIS 4. BL 
1 — vp, ay Ak l— »(aga&— l)’ 1— ak ae RARE vy (ak az + 1) 
D. Oe LE ee ee TERN 1 
i 1— vz (ai ar + 1) À ara 1— v&(a& az — 1) wg iun DEM E T 
| (2.21) 
_ With eqs. (2.21) we obtain 
lanad UC UN ee 
hence 
| = ^ PN ^ ^ ^ à A h fed 
| = [Pr B] =0, [fe à] = 79x (2.23) 
| The subsidiary condition are now 
(i, 4990, SR come (2.24) 
| or 
| | Er a-p + a& (1 — vk ak a | ð= 0, Ik] S ED (2.24b) 
\ Tl—va rar 
and, finally, A 
OU" 1 la, 4- (1— vx a^ a. y) a, (0 — 9x ax a,)| 9 = 0, PES Ge 
We put = Fe 


OH H 9, (iT) 
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We cut off the infinite sum in (2.25c) putting 
Dp = ar Sansa (2.25e) 
Nk . 
wherein n, is a natural number. Substituting (2.25e) yields 


l 
(az s) t]*s 


nk 
Dy = D = o. (0,0) (ma? Y. (HI) — r (az at)! Bo (2.26) 
ix k 1 


s=0 


or, on normalization 


n2^k 1 
"k — 2h te OS Ww tcI E uoce “or ys l 2, 
^ ht i > CD an Ve 070 


Thus it has been proved that the original formulation of the Bohm — Pines theory 
can be approximated to any degree of accuracy, provided the oscillator operators 
are replaced by their nearly-Hermitian analogues. The theory proposed differs from 
that given in (19582), since the equivalence with Schrödinger’ equation is now direct 
and not on average. This is made possible by (2.23). 


3. The partition function 


In order to derive the correct statistical function, we apply the identity, Dyson 


(1956), 


Bm-1 


60 b B A : d 
eB AL +Å >: 2; D” f df f df, + f dp, e FA, . eh P, ehr s bH, e` fH x 
m= 0 0 0 


x … x en H,e-fmAı, TT À] +0. (3.1) 


The identity is proved in Appendix 
We choose the extended. Hamiltonian and the subsidiary conditions according 
to the author's paper (19583): 


^ ^ 1 ^ * 
Hexe = H + 4 p» ho, (aja, + a,aï + a,a. + aza) — 


5 
la] <a, 


1 == 
= hi ho; (a_; + ai) a 4 5 Dos ho,, . (3.2) 
es = A 
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where the x, are renormalizing parameters 


Oscar, L (3.4) 


and 


In the foregoing relations À is the reciprocal lattice vector, œ, — the frequency of the 
collective oscillations, 


An ne? 1 1 Z l 
2 —— = ——J — 
(03 = Wp = Are — eam TE (3.5) 
n — the number of electrons, Z — the number of valence electrons per atom, m 


and M — the masses of electron and ion-core respectively, 9, — the Fourier transform 


of charge density and H the energy operator of Schrödinger’s equation 


n N 
pe Pi pz 2n 030-1 — (1 + Z) ne? 
die OX doe i (3.6) 
=. = KO 


The total charge of the system is assumed to be zero. 


We have further 


Dam ee Dy ar Tasche | (3.7) 
Al<a, 
4,70 
with 
et Se Conn Tutor (3.8) 


Applying the subsidiary conditions (3.3) we now transform certain non-diagonal 
terms and the interaction member of the Hamiltonian (3,2). This yields 


a, = — x,2; 9, a, a. ., D = — x, a, a; 9, (3.92) 
T * * n * 7 

a*., D = — — a, D, aiat, = — — aï ai Pa. (3.9b) 
Xa XA 


The applied transtormation is permissible, when the difference 1—x, is suf- 
ficiently small, since the Hamiltonian (3.2) and the operator of subsidiary conditions 
(3.3) do not commute. 

With the relations (3.9a), (3.9b), (3.9c) and replacing w, by ®, we have 


ER SET BR hern c i) 
Rea = À + Mn y too —ad stent Duy uk mes 


= ii (<a 
ji lal&Ae. ; le A, 2,50 (3.10) 


1 
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LA om 


OUU RP eee 
wherein 4 
2 W 
Q, = o, fie LE, c1) 38 
= 
It is convenient to divide the Hamiltonian (3.10) into three parts a 
Hu. = Hy, + Hy + C, | (3.122) 2) 
with 
RISA 
He Dh; ho, (1 — = “0 i (3.12c) 
_ HA, : | 
= ho, XA (l = 2 rans 
IA, 
220 


We can now proceed to calculate the partition function, By (3.2) the trace 


of exp (Bit a) (where 8 = 1/kT, k — Boltzmann’s constant, T — the absolute 
temperature) is computed by averaging the Enc y over the plasmon part and | 
taking the spur with help of the functions y (... Tj ...;... t, ...). Thereby 


Z= ae tas, 813) 
where the mean I denotes; sy zn ! | 
ze) dern 

ae d ea ee Her Are 


mar a. 
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We have then 
Bm-ı 


een eo ED = ROT (Ur Jan f 48, «e (B,) x 
m=0 0 


E x n B 
x Hy (D... Ely (By) > = ettet f dB, ce BS) + 


Bi 


B Bs "T A 
+ f aß, | dpa e H, (Bx) Ha (Bd) > — f df, J df f Jfy «e A, (61) He (Pa) X 
0 0 


Ba 


p B Bs 
x Hy (Bs) > + f df: | dB: J dhs | aß, < e PP Pr; (By) Hr (Bs) He (Bs) X 
0 0 


x Ha (B) > — +} (3.16) 

wherein 
Ê, (B) = «PH, e7. (3.17) 
The relation 
o, ay = 0 

yields 

D, = "m d» e BA - (3.182) 
hence 


> PAT = EE 


CNE, 
x H, (B3) H, (Ba) -- "Ha (B,) H e 9, 02720 Dy, (3.18b) 
Ir 
with 
a% (0) = e° Maj retha (3.18c) 


The following relations are readily proved: 


[aa (Bj). a, (B,)] = 077 AC (3.19a) 
and, quite generally, 


min (I, s) 
a, (8) la; (BY = » () ( ) ZI elr Por? [a5 (By) 7" a5 (B;)- (3-19) 
^ r=0 
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The first term of the expansion (3.16) is now 


<e Ph = 
=e? IT (x) d, expl- % x, 4,(8)a_,(8)] expl 2j x 4% (0) a~; (0)] 6, — 
<A, CEZ LESA 
4,>0 A,>0 0 
= PH IT (1— x3) dt exp [— x, a, (8) a. (8) - exp [— < a} (0) a! (O)] 4s = 
R<, 
4,>0 


er) Ay La D$ a} (B) a... (B) [a5 OI [a*., (0) D, = 
[al<a, T=0 


co s oA 2 
=ef IT aay € ae I! s!6,,e7PSh 22 = eb IT a 


75 2 B= BE One 
jal & 4e 1=0 Is! [A] & 4e l—xe PRO 
A4z>0 s=0 Az>0 
(3.20) 
The second term vanishes; 
eg H: (04) > = 0, (3.21) 


as all matrix elements with odd products of H, (B; vanish. This results from the 
fact that the plasmon part of the function (3.7) has the form of even products of cre- 
ation and desctruction boson operators. Quite similarly the even products of a; vanish, 
provided the sum of index vectors À is non-zero. 

In order to compute the third term in eq. (3.16) we note that 


H, 2 (Bs) Ha (Bs) = > 2 h,h,(1—x,) (1—x,) (hw c)? [a3 (B1) “(BD +a", (Ay) 8 s-4 (81)] X 


payee 
X [a (Bs) Co (B2) + a (B2) 0-5 (Ba) = Dy Hal — a)? (ho)? [a; (By) a* (ba) x 
Äl<A, 
A720 


X 04 (B1) 0-2 (Be) + a, (b1) a (Ba) 0-4 (B1) 0, (83)] + ... (3.222) 
The terms marked by dots vanish on averaging; hence 
= 1— x 
AG Ba) = AU s BAA (TO x 


His 
[7l <%e Az>0 


x die PP TAT SE Ar 559 ds t5 (B,) B-0(Bs) + 8-0 (8) 8, (B,)] = 


1— x xg e(&Bi—28)80,2 
= —¢ -PH (howp)? p» hi (1 = Xv) (1— xo)? | JR a à Ko e eh PD: by 
[6| € 4c DE 1— x eho ae): 


oz>0 


X [s (81) -o (B2) + 6-5 (b1) ĉo (Ba)]- (3.22b 


1 
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Finally, we have 


ds K l—x hi xs (À — xo)? 
BH, ne. AS ARR Acid c ; 
>= PT rau. hot ST reine * 


Hx A ió 
Az>0 lel ^e 
x E(B: + B2—2B)A2cI2, 6, (By) O_o (Bs). (3.23) 


The higher expansion terms, which are somewhat more involved. are of similar 
structure, and exact investigation proves that they are not necessary for the present 
computation. 

We have now 


2 hz o Xo ELI 
<e bH — nle —6C II He er en a — (hog)? yia Lo X d = 
e e 


ia 
| 3 Y jo] <Ac 


ß g 
x eine [ dB IE „erhal? 6, (By) ebifi2cl2 6... (By) + Jr (3.24) 
0 


0 


Let us investigate the form assumed by the equation (3.24) for the Bohm — Pines 
(1953b) subsidiary conditions. We have: 


x; l — 0» Qe <ul (3.25) 
ho (0 xa) ho, ĝa (1 — 01) | 
= : - x 3.26 
c 4 Qo 1 ve > 2 2—ô, ERZ 
AER HEL 
A770 1,70 


For sufficiently small ô}, 
ho, 5 (1— 6, F) 
re py PFARRER SEXTO COE (3.27) 


The quantity (3.11) in this case assumes the form 
14 X ee 4 c 
2 =o, (1- £3) - o [17 za 2 


li l 0) 
of tosara+-2+ 20228 (3.28) 


The expression 1-21, 2 e-Ph23 transforms to 


Eye Sar 2/1 — 0) 
1 — xt = 1 (l Que i P e Se z l4 pho, 64/4. 


PEN MAT SS ee = 
Er E em 2 Oat Bho/2) Ò 
3 Xi € 2 1 a) (3.29) 
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3 40, 2 r 
CRE || df, fa, e^ 73^ 0, (By) oe 2 057 (Bo) i df, f dB 0, (5) 0_ (Be) 


0 0 


(3.28e) 
and 
Run 62 (1 — à 
ptc = BL) ur BE SIE i (3.31) 
1 — x2 eho 26, — (1+ Bhw,/2) 62 2 


We have now 


^ ANA -2 
« ei) un > = oP IT ef? E: 4 


lal<a, 


40 


x (L + Bho, 64/2) x 


B Bh | 
x (ho) 27 Med | dpa f as êo (Bs) êo Ba) +. = 
9| ^c 0 : , 


B Log 
= eee (1—f howp 44/4) (1+ B herp 64/4) [1 (ho) 20,12: k df, / Bx ĉo (b1) x 
Apo 


Xê- «(B --...] = ui h- (hein)? » hi9, f jaf dfs àc (b1) +. ] 


<a, 9 


For sufficiently small do one gets 


- 


qu ur D | l 3 


ra 


Le 


25: oM 


ee gs 
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We divide the Hamiltonian (3.34a) a gain into three parst (o; = ®,) 


We HH EG (3.35) 
H, ZH + has A Dp ho, (a; + a. 4) (a, + any, (3.35b) 
<a, 
= ye h, ho, 0, (a, + 25). (3.35c) 
<a, 
Gers UE (3.35d) 
<a, 


Since the plasmon part of the state vector is now the vacuum vector ®, the 
averaging is very easy one giving 


EUER S EO AU MPH AIRE TU ofl, 


+00 B Bm; ^ 
x D c [ 4B J df, Hy (B)... Ha By) Po (3.36) 
m=0 ò 
wherein 
Bi, (8) — Ah pehtm D h hop âa O) Ca tad. 630 
Â<A, 


We have further 


B y 
e au = e KC 9H) [S e Ph, 2 hy ho, | dpi ôa (£1) Dy e P^ x (a.., + aj) Do + 
A 0 


Al<a, 


ar 22 h, ho (ho)? f df, f df; 0, (P1) Go (b2) D e-P^(a. , + a) (a, + a a) Bo t ..] 


«A, 
jal<a, (3. 38) 
Let us compute the first term in the expansion (3.38). 
+ pid = II Die —Bho,/2° - (a. 42 laatat 2 D y Che y 
| y " «A re zer 
Az>0 az>0 
| ^ (bhos?) Gi or 
x Di (a-a + a) (Ca + at) Do = H T T Dé at; (ai Do = 
fil<a, s=0 
Az>0 
Supe ho io qp Lage (3.393) 
== II ! Ss: <a, 1 + Bhop/2 


AA RES 
Rica 170 Az>0 
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GA 
Ke] 
O 


Bho,/2 < 1. (3.39b) 
We have used the equality 3.19b) (8; — f; — 0) and the relations 
AUR uar = 


The condition (3.39b) is not satisfied in low temperatures. The way to overcome 
this difficulty has been given by Dr A. Pawlikowski. The following computation 
belongs to him. 

Let us write the destruction and creation boson operators as follows: 


1 ee 1 S 
a= Vi (a; + iai), af = yi (ucc ia), (3.402) 
Qc Se TM (3.40b) 


where a7, ai and a/*, ai* are the real and imaginary parts of destruction and crea- 


tion operators. 


Since, according to (3.40a, b), 


(aa le eos ei (3.40c 
[aa Drama] = EE (3.40d) 
the Hamiltonian part (3.35b) may be written as 
f _h 
a B» (a, +a +) (a_, Se at) = P CA D =) (ar Je art)? (ai ST ait)2, (3.41) 
u |a, s BE 


When we denote the plasmon vacuum vector as |0>, we obtain instead of (3.392) 
the following expression: 


« oje^5 o. = = ID <0le ara Mara go. — 
A<2, 
RC 
zv IT SRE [0 Ofer a? 0>, (3.42) 
il< a 
2,>0 
with 
hw 
oki (3.43) 


Each of the square terms in (3.42) may be brought into the linear form by means 
of the well known formula 
+00 
d l Ber 
eu — —_ [dxe-x-uVe xu, 


IT 
—oo 
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We have then, 
+00 


+ il NAA TEE 
< 0|e-7(5 * 410 > = = Lx « 0| eti Vox aj + 232 |0>. (3.44) 
7 


On the other hand two boson operators satisfy the relation 


+ 2 c 5 
ea(a+a jy e et le. ewat «eX. (3.45) 


In connection with the last formula we obtain finally 


+00 
< |0 etaa 0> = < Ojeli tart |0 > = a Í dx: e 4-299 x 
TT 


— 09 


+ oo 
7 2ix Vo a5 5 QixVo GA = ui. —(1-4-2e)x* — 1 
x =< Ole are 110> ys fe (L3 299 (3.46) 


Returning to the equation (3.42) we have 


1. P 1 
« 0| e- Ph. 0>= II =) ll IT TOP (3.47) 
Asa, \W1+ 20 «a, LT 20 i cape 1+Bho,/2 
Az>0 Az>0 nd 


We have then obtained the same result as in (3.393), but without the supplemen- 
tary condition 


pho, 
ITE l^ 
2 m 


In the expansion (3.38) vanish all odd terms. The second non-vanishing term 
has the form j 


h B Bi s 
y, (hot [ dB, [45.0.69 dx Ba) Poe À (2-2 + a9) (ea 9o = 
0 0 i 


Als, 


= II > itae, f ap f dada “a (Ba) Po e Lite De à I 
<a, MIA 
X270 


- (a-, + a5) (a, + 83) Po = 


=2 > (ho)? | df, j IRR do AY CRGA Li AE NE EN 
fl <a, 
Az>0 


# 
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qoot: l 
x (Gy, A a?) (a, me 2) TI e` Bhogl? - (a gt %) (ag+a— 9) D IE T+ fho? X 
[e| & ^, lel «4, 
P 0z>0 
oz>0 
A B Bi 
h? (ho)? ) li E ^ 
x ——— | d d, ô— (Bo). (3.48) 
Yon) sf 46 6000: 
l< 0 0 
From (3.38) and (3.48) we obtain 
_ppext (c+) IT 1 14 h(ho)? 
Dux prd Edu NS cer T v7; 
Spec ; [l<e 
B Bi 
X f a8 | ae ôa (B3) 0— > (B) + E (3.49) 
0 0 
The partition function assumes finally the form 
IUE ^ 1 
Z= Tr < e PHM > = e-8C Tp e-BH. [J —— | d 
pese PER 1+ Bpho,/2 | 
oz>0 


TT 
hi (hep)? A en } 
=> T+ Boy dfi ji df, ĉa (P1) 0-1 (BD + | (3.50) 
Al & 4e 0 0 


We see that Kazanawa's theory (1957) leads to the false statistical function. 
. Should we define the renormalization parameters as 

l kT 
=1-— = 1- ——, 3.30 

= pho, ho, 
this would suggest that the formulation of collective description given by Bohm 
and Pines (1953b) is valid for sufficiently low temperatures, whereas the approach 
of Kazanawa (1957) extends into range of temperatures wherein heat phenomena 
(electron collisions) play an important part; however in the latter circumstances the 

formalism of the theory becomes non-linear, Bohm and Pines (19522). 


4. Conclusion 


In section 2 a new renormalizing procedure is formulated. The variant of col- 
lective description proposed by the author is nearly-Hermitan and approximates 
the Bohm and Pines theory to any desired degree of exactness. Equivalence with 
the many electron Schrädinger equation is direct. | 

In section 3 an exact statistical function is derived. In connection with the comp- 
utations it was found that the Bohm — Pines subsidiary conditions, improved by 
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introducing a certain renormalizing term, yield the exact partition function, whereas 
the Kanazawa’s supplementary conditions lead to false results. The possible inter- 
pretation of the renormalizing parameters as functions of absolute temperature 
is suggested. 

The author wishes to express his best thanks to Professor R. Ingarden for his 
inspiration to enter upon the present investigation and his scientific coworker Dr A. 
Pawlikowski for the computation concerning the partition function in the case of 
Kanazawa's subsidiary conditions. 

The author is also greatly indebted to Professor Dr S. Szczeniowski for his valu- 
able help during the preparation of the present paper. 

I thank also Mr R. Ferchmin for explaining some interpretational aspects of 
Dyson paper (1956). 


Appendix 
We prove the validity of (3.1). We have 


s—2 s—2—l, 


(H, + H,) = Hy + z OS A: RE 2. > Hi À, Hh H, H:-3-0:09 + 
=0 l= l= 
—k s—k— s-k-(h+...+k-) , RR ^ ACER A 
+5 = V > Hh A, Hh H, ... Hit H, HA 0-9 E uuu + Hi 
1=0 4-0 k=0 3 


(A.1) 
Further 


e&t. B) — er? ip E fnis er? (lp. E p SA Bb A mt a 


120 
ARLAN EN 
= HY A, Hy He H-2- +1) FP oo AF ~ 


MEL i e TEN Se MORES ane 
re s ; >> e Y BL DS Hy Hi gt a 


h=0 h-1 Ik=0 


E d (— 1)s+1 pst! z a p, fe^ 
+..= may. men | Ahr 


s—lh 
eine ih B.E À, ETT E uu 
1 2 Í 2 eee 
+5 63) cea 2, 

s+k As+k CALE a ^l r? th D She py (ht... tlk) 
jones B y Eh = Abo, Bb B. ÀF À, È EA 
GB h=0 i,=0 =0 (A2) 
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On the other hana, 
en +00 B Bi Bm-ı paie: ze 
e KA Bo — y (—1)" f dfi [ 46... f dp,, e^ Ps Hy e(bs— PF H, X 
n0 0 0 0 
— " cal ee dd Sp) d Pi 1 
x em Bm- f H, g- m Hi = g- Hi — f dp, PDG Hye Ph + f df, ji df, ei Hi x 
0 0 0 


B Bk-1 mA A 
x Aer OH, el. +) Í ih [E dp,e® PP. Her Pr Het... 
0 0 


(A.3) 


The first terms in expansions (A.2) and (A.3) are the identical. Let us compute the 
second term in (A.3). For this purpose the following integral is required: 


Bi 
A $ 
il die Min CN ar NO ESS 
0 
Then 
f E V vcn ^ 
= [ee e(dı-MHı H, efi = M y si fers Hi H. Hi E 
ò l=0 s=0 Ü i 
+ oo + oo 1289 z 
T AVET N (—1)"+1 AE vn. 
= NS H H. H, = n+l E nr 
u He IRD ie e (ni^ 2 mA 
$ n= s=0 
A (— yn TA quias: 4 
= Sr YaBB. (A.5) 
s=0 |. 4-0 


The term derived is identical with the second term in expansion- (A.2). 
Quite generally, 


+ oo + oo s—(h+...+lk-1) 


+00 
; (— 1)s tk ps tk nl at wh, 13 ^ SES 
Hr H, H? Ĥ, is Hit H, H (h+...+Ik) < 


s=0 Gen h-0 1,=0 Ik=0 
B Bi Bk-1 | 
zer J df, J df, ... J dB cP fi, et Er, ... P mft FT, Pr, 
(A.6) 
The equality is proved step by step. According to (A.4) and (A.5), 
Bei + © n gntl |" 
[ dp, em hi-fi. Hy eom = CDD T YRA (Aa) 


0 n=0 p=0 
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Furthermore, 


Bk-2 Bk-1 
f dB, 1 f dp, e(Pr-ı-Pr-2)Hı H, eBk—Br-1) H, H, el s 
0 0 


+ oo Bk-2 


| ROME RENT 
— SM f dpr- mo 1 e(Br-1—Pr-2) À y. IR He AR? = 


=0 - 
n : 0 


Bk- 
ry 22 SH dßr-ı (B1 — Br—2)! Bani Hi A, Y at fe 


p=0 


= EC ees Or in m À, HP À, m 
2222. entr m À, Amt À, Hi 2 


We put n+s = l, obtaining 


Bk-2 Bk-1 


f df, 4 f dB, e Pr-ı-p-DAı FT, Pe Pr- nét, jf, ee — 


+ oo I 1—s (— 1)! 
a Y Erz BSH IT I DS en . (Ab) 
1=0 s=0 p=0 ; ^ 
The third integration yields 
Pr-3 Br-2 Bk-1 
x f df, . A f dB. = J df, e &k-2— Br-3)H, H, efe 1—8x-2)H, H, e &&— Br-DH, H, gil. — - 
0 


Tło $ s-h sch- quaro 1y 
EX > 2 » EE 8:43 Ht À, Hb Bert, (AZTO) 
s=0 h-0 h-0 1-0 
The: k-th step gives the identity (A. 6). 
Let us now investigate the form of the fourth term in (3. 16) (the third term van- 
ishes). For this purpose we select from H, (Bx) H, (Ba) H, (Ba) H, (B4) the terms, 


which do not vanish on averaging. Let us write them symbolically as follows: 


23-(2+2+2)(2+2+2)- 


ES x-—— x. 
=e e E60 


M 
ar 


on x~te 
= ERS > + » 2 te (A.8) 
e Ne e na E jme, 
E dx9 x+Le 


J. Sledzik 


02 
Hence, according to (3.12c) 
Bi, (B1) Ê (Bx) Ba (Bs) Fly (Ba) = (ho) I 9; 0— —x,)? (1—x,)? [a4 (By) a5 (Ba) ax (Ba) X 


[|< de 
ae 
x a, (Ba) 0,(Bı) AEA, 0, (Ps) e, (Ba) + a(b) a” Xf) a (Ba) a. (Ba) 04(P1) 0 (B3) 0,(Bs) X 
x6-.(B4) + 23 D hy lg hy h (1—21) (1—20) A) (17, ) a (Bx) (B2) as (83) a(b) X 
PLACES 
Joi <A lel< 4 


4g xe prm 


x 68) 63) Ox(Bs) Go(Ba)} + ++ = (hop À > hla) [a5 (83) aBa) (83) X 


[4| & 4e 
X a. XB4) 0,(Bı) 0 (Be) 0 ,(Bs) 0-1(84) -Fax(f1) a* (B3)a;(B3)a. (B 4) 04(8:)0- (82) 0,(83) A 
x e (Ba) + a; (Bı) vu (Be) M (Bs) a, (B4) 0; (B4) 0- a (B9) 6 a (Bs) 0 (84)] + 
AE a hi hy (1—x,)? (=)? [a; (£1) a ,(Bs) a,(ß3) a, (B4) e£) 0-1(B2) 0,(B3) 0-,(Ba) Pr 


jÄl< àe 
Ixl<Ae 
i++r 


+ a, (£3) a, (Be) a™_, (Bs) at (B4) 0 (B4) 0, (Ba) o- (Bs) e. (Ba) ar 
+ a3 (B1) a5 (Bs) a*s (Bs) a*a (Ba) 6a (Bx) 6x (Bo) ô- (Ba) 0-,(B)D + (AI) 


The procedure of averaging is an easy one; we give the final result only: 


f ae, f dj, j df, f dj, < ef Ê, (B) Hy (Bs) Ha (Ba) He (B > = 
0 
^ Tto ho xs — xc f f a 
= qs — err eho > fe rue) j df, | df, L df, | df, x 
LES 


x oet Aug 40) 0012 [es (B) à, (Ba) 0. (Ba) 8-0 (B4) + Go (B1) 8-5 (Ba) ĉo (Bs) 8-0 (Ba) + 


| + c (£1) 0. (Ba) 0. (Bs) 6« (B4)] + (hoy). Aen na, 
ash 
dx 


B Bi Bs Bs 
x f dB, [dB f ds f ABa [etA hoe, utet t 2PRONA- OB) 0. (Bs) 0,89 0- (0) 
0 0.2 —0 0 A 
+ ePıtBs-2B)hQe]2 . ePr+Br-28) ROM G 2, (B3) 6,08 2) 0-, (Ba) 6. (Ba) + eBi+B—28)h0012,, 


Ken AO 
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By a procedure identical to the one of the text, the term (A.10) is proved: in 
the case of the Bohm — Pines subsidiary conditions to contribute but negligibly. 
The higher expansion terms in (3.16) behave quite similarly. - 
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EINFLUSS VON ÜBERNACHSTEN NACHBARN AUF DIE ÜBER- 
STRUKTURELLE ANISOTROPIE DER SPINWELLENDISPERSION 


Von HENRYK COFTA 


Institut für Theoretische Physik, A. Mickiewicz Universität, Poznañ und Abteilung 
für Ferromagnetismus des Physikalischen Instituts der Polnischen Akademie der 
Wissenschaften, Poznan 


In einer früheren Note (Cofta 1957) haben wir auf das Problem der Abhängigkeit 
der Spinwellenausbreitung von der Überstruktur des antiferromagnetischen Raum- 
gitters aufmerksam gemacht. An Hand der angenüherten Rechnungen haben wir dort 
gezeigt, dass für manche antiferromagnetische Überstrukturen eine ausgesprochene 
Anisotropie der Spinwellendispersion auftreten muss, hervorgerufen ausschliesslich 
durch die Überstruktur. Hier wollen wir einige Resultate der genauen Rechnungen 
vorstellen. 

Wir kónnen uns auf die Formel (25) von Cofta (1959) stützen. Ihre Anwendung 
auf die einzelnen Überstrukturen der kubischen antiferromagnetischen Gittern 
bestätigt unsere Vermutung, dass die überstrukturelle Anisotropie der Spinwellen- 
dispersion in allen Fällen auftritt, wo die parallele ferromagnetische Ebenen - mit 
abwechselnd gerichteten Spins vorhanden sind. Der einzige Ausnahmefall betrifft „ge- 
wöhnliche“ oder „natürliche“ Spinordnung (d.h. alle nächste Nachbarn haben ent- 
gegengesetzte Spins). Die Rechnungen zeigen, dass das effektive Feld H4 der magneti- 
schen Anisotropie fast keine Rolle in er unteruschten Erscheinung spielt. 

Die Berücksichtigung von übernächsten Nachbarn führt zu sehr interessanten 
Folgerungen. Da die Wechselwirkungen zwischen den übernüchsten Nachbarn im 
allgemeinen einen anderen Charakter haben als die zwischen den nächsten, hängen 
die Anisotropieparametern von dem Verhältniss n = JılJa ab, wo J, das Austausch- 
integral zwischen nächsten Nachbarn und J, — das zwischen den übernächsten be- 
deutet. 

Diese letzte Frage wollen wir an einem ziemlich einfachen Beispiel verfolgen, 
nümlich an dem FC-Gitter mit (001) ferromagnetischen Ebenen. Aus der obengenann- 
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ten Formel erhalten wir in diesem Falle für die Kreisfrequenz der langen Spinwelle 


(ak 1): | 
w = |y| (C + 2a? [H 4 (H, + 2H;) — 32H, Hj] (42 + k?) + 4a? [H4 Hs + 
+ 16 H, (Hy, — Hj] #3} () 


Das Symbol y bedeutet hier den magnetomechanischen Faktor, a ist die Gitter- 
konstante und £,, k,, k, sind die Komponenten des Wellenvektors. Ausserdem sind 


in (1) die Abkürzungen 
H, = S Jjyh und  H,-— S* J,[yh 


sowie C = H,(H,—32H,) eingeführt. Da das Feld H, bekanntlich sehr schwach 
gegen das ,Austauschfeld^ H, oder H, ist, kónnen wir H, vernachlässigen und 
schreiben: 


o = 8|yl a (H, (H, — Hy) k? — H, H, (k? + k^ (2) 


Für die Wechselwirkung zwischen nächsten Nachbarn (J, = 0) erhalten wir eine aus- 
geprägte Anisotropie: 


o = |yl {C’ + 2a? H, [H4 [02 + k? (— 32H, E] (32) 


die die Richtung [001] senkrecht zur ferromagnetischen Ebenen auszeichnet. Wenn 
wir dagegen nur die übernächste Nachbarn berücksichtigen (J, = 0), dann wird die 
Dispersionformel isotrop: 


o = |y| {C" + 4a? H, A, ky (4) 


was sich leicht auf Grund der geometrisch-physikalischen Betrachtungen erklären 
lisst (Cofta 1960a zund 1960b). Beide Effekte sind demnach sehr verschieden. Die 
isoenergetische Fläche E = const stellt für die wahrscheinlichste Werte (Cofta 19602) 
1j «0 ein Rotationsellipsoid dar. Das Achsenverhältniss a:c hängt in der Näherung (2) 
nur vom n ab: 


= -y1-m à (5) 


Die Abhängigkeit (5) schafft die Möglichkeit einer unmittelbaren Bestimmung von 7j 
durch die Messung des Parameters a:c der überstrukturellen Anisotropie. Da wir 
bis jetzt fast nichts über das Verhältniss 7 in einzelnen antiferromagnetischen Stoffen 


wissen, wäre eine solche Auskunft sehr wertvoll. Die móglichen Messmethoden werden. 
spüter bearbeitet. 


Die analogischen Untersuchungen lassen sich für alle antiferromagnetische Über- 
strukturen mit ferromgnetischen Ebenen durchführen (Cofta 19604). Als eine 


' ^ 
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sehr nützliche Hilfsmethode erwies sich dabei die Einführung einer geeigneter Klassifi- 
kation von Überstrukturen (Cofta 1960b). 


Die Einzelheiten der genauen Rechnungen sowie volle Diskussion werden in 
Kürze in Acta phys. Polon. publiziert. 
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NON-LINEAR EQUATIONS FOR SPINOR FIELDS 
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A plausible form of a non-linear spinor equation, similar to the equation of 
Heisenberg, is derived by using the principle of correspondence with general relativity. 


Several authors suggested that the exact equations for spinor fields should be 
non-linear. In particular, Heisenberg(!) investigated two versions of such equations, 
based on the Lagrangians 


E je - piz - 
Pu On V t. (py)? or wvy,9,y tz; (v 7,75 p)’. (1) 


Such assumptions are usually believed to be merely examples of the simplest kinds 
of non-linearity, whereas the exact form of non-linear terms is not known and may 
finally turn out to be rather complicated. However, criteria for discriminating a parti- 
cular type of non-linearity are not yet known. The present note is devoted to the 
problem of determining a plausible type of non-linearity by advocating the principle 
of correspondence with general relativity. à 
At the present stage of development of theoretical physics gravitation plays 
a distinguished role. Being completely amalgamated with geometry, the theory of 
gravitation is free of arbitrary and simplifying assumptions to a degree surpassing 
by far the o'l:er theories of fields. And just the equations of gravitation are non-linear 
and the form of their non-linearity is well determined (except for the disputable 
question of existence of the cosmical term). This fact constitutes a strong argument 
in favour of the opinion that — similarly as in the theory of gravitation — linear 
equations constitute only the first approximation valid in the limit of weak fields. 
At the same time it constitutes a hint as regards the problem of determination of the 


_ proper form of non-linear field equations since there may exist an analogy between 


the type of non-linearity encountered in the equations of gravitation and in the equa- - 


tions of other fields. : | ` 
Since (covariant) lagrangians for different fields have to be added to the lagrangian 
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of gravitation, there should exist some relationship, some analogy between these 
lagrangians and the curvature, otherwise their union into a common lagrangian could 
be regarded as an artificial *cementing of quite different pieces". 

In the unitary theory of a gravito-electromagnetic field (2) there exists indeed 
a close relationship between the lagrangian of gravitation and that of the electromagnetic 
field: both appear to be parts of a common lagrangian being the curvature of a six- 
dimensional space. However, one must not hope that by a further extension of the num- 
ber of dimensions of the underlying manifold it will be possible to include all the 
fields playing a role in Nature. In particular, a lagrangian with the simple geometrical 
meaning of curvature cannot account automatically for the existence of spinor fields. 
Therefore the relationship between spinor field and gravitational field lagrangians 
must be of a different character from the relationship between the electromagnetic 
field and gravitational field lagrangians. 

Since the curvature is expressible in terms of Christoffel symbols and their 
first derivatives, the analogy between the lagrangian of spinor field and the curvature 
may consist in the fact that there exist some analogues of these symbols (let us call 
them I7, [y]) and their first derivatives (let us call them 9,17, [y]) and that the lagrangian 
for the spinor field y is the same function of these analogues as the curvature is of 
Christoffel symbols, 


L(y) = Cg" (9, I$, fy] — 9, T5, [v] + Ti (vl TR. lv] — Dos [yl Tr, vl). ©) 


The fact that the lagrangian for a spinor field is to be obtained (up to the factor C) 
from the expression for the curvature by mere substitutions of the Christoffel symbols 
by their analogues, explains their relationship and the necessity of fusing them into 
a common lagrangian describing gravitation together with its sources in a more or. 
less uniform way. 
Since particles with half-integral spin admit two concepts of velocity given by 

the operators 

> Re TS fa 

a= iyty) Lande. Ga oe (3) 


(to be generalized in a straightforward way for four dimensions), the expressions 


v m vo 
Yu — wy EU Y! Va Sua + Yu Bie — Yo Siw) (4) 


constitute some analogues of the Christorrel symbols. More exactly, the analogues 
of the Christoffel symbols we are looking for should be proportional to the densities 
of the above expressions 


I [y] = ky yip and OT, [yl = Aya av | (5) 


where ô, operates upon everything to the right, and k is again a factor of proportio- 
nality. 
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In contradistinction to the usual Christoffel symbols and their derivatives, their 
analogues (5) do not vanish even in the limit of flat space and in inertial coordinate 
systems. Thus, for reasons of simplicity, let us neglect the curvature of space-time 
and consider the expressions (5) in an inertial coordinate system, introducing an 
imaginary time coordinate and the usual (hermitian) Dirac matrices. Under these 
assumptions (5) assume the form 


Dv] = x (Pa 9, + Yu das — 9,94) V (6) 
2, D [v] =% v (y; Su» Tr Vu Oxy br On) 9, yp. (6°) 


where y = y" y‘ and x is again a proportionality factor. Introducing (6) into (2) we 
find the plausible non-linear form of the spinor field lagrangian 


L (y) = 6Cx [((vy,9, v) + x (oy, v) (ry, v)] (7) 


where 6Cx is the coupling constant of gravitation and k is another constant with 
dimension [cm?[ in natural units h = c = l 

Thus, we got a lagrangian resembling closely that of Heisenberg except for the 
factor y® in the non-linear term. 

It is easily seen that the possibility of supplementing this lagrangian with a term 
of the form myy corresponds exactly to the possibility of supplementing the curvature- 
lagrangian with a cosmical term. 

"The essential step in our procedure :was the replacement of the derivatives 9, 
in the Christoffel symbols by — my? y, according to the two concepts of velocity for 
spinor particles, and a subsequent transition to the densities y" > y" yj, V. Such 
a procedure can be understood intuitively in the following way: there are two aspects 
of the motion of the spinning particles. One of them, connected with the Zitter- 
bewegung, does not agree with the classical requirement of the geodesic motion of 
a free particle. But any departure from uniform motion of a free particle can be 
understood as an encroachment upon geometry. The expressions (5) determine a sort 
of “affine connection” accounting for the “geomtery” where Zitterbewegung replaces 


the uniform motion along a straight line. 


(1) Heisenberg, W., Rev. Modern Phys., 29 269 (1957) 
(2) Rayski, J., Nuclear Physics, 16, (1960) 


is 
> 
. S 
* 3 ~ 
- LT 
dd Line 
e 
* om A [oou ed d 


rs De A pH 


E areas, ripped 


"e CL m adi Tui 


denis qu den; 


+ 


RUP 


x 
k r 
te 4 x 
: Eb 
2 - 
e 7 zw #3 
n Y. E 
= ` „El: wid 
: 5 ke 
Di ~ » 
ww , e 
g » k 
$n 
H a. 
S 
+ wr 
4% Y 
Sr} T3 u IX ^ à ` 


M. pe "yg rn | GN. io. 


Pte nl "RES ar =e = dt E 
Ann yon sid n 


| 
| » 


Volumen XIX — Fasciculus 3 


Rytel M. Sur une méthode d'évaluation de l’anharmonicité des 
molécules polyatomiques . . 273 
Bogdanowicz J. Danysz M., PUDOR O Mare 
quit E, Skrzypczak E, Wroblewski A., and Za- 
krzewski J., Determination of the Mass of the A? Hyperon . 277 
Fechner B, On the Statistics of Spin Waves by the Bethe Method. 289 
Kopeé Zb. Effective Mass Method in the Case of Non-Quadratic 
Dispersion Formula . . . Le 200 
Ostrowski K, Rectification of Photoelectrie Spectron 25 3:219 
O'Connor D, Sosnowski J., Measurement of the Slow Ne- 
utron Spectrum of a Neutron Beam from the WWRS Reactor by 
Means of a Crystal Neutron Spectrometer . . . . . . . . . 329 
Ostrowski J., Collision Recombination in TS. . . . . 839 
Klimowski y. and Pietrzak J., Dielectric Permittiviiy Varia- 
tions in BaTiO Single Se and Ceramics Resulting From 


Hydrostatic Pressure . . . . 869 
Sledzik J., Some Remarks on Collective Description of Meinen 


Interactions. The? Partition Function. e e a n e 0$ 383 
Letters to the Editor 
Cofta H, Einfluss von übernächsten Nachbarn auf die überstruk- 


turelle Anisotropie der Spinwellendispersion . . . . . . . . . 405 
Rayski J., Non-Linear Equations for Spinor Fields . . . . . . 409 


